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FOREWORD 

The SHABERTH computer program was originally developed by Kell- 
strom [11 under U.S. Army Contract UAAD05-73-C-00U sponsored by 
the Ballistic Research Laboratory (BRL) , Aberdeen 'Iproving Ground, 
Maryland, to simulate the thermo- mechanical performance of load 
support systems consisting of a shaft supported by up to 5 rolling- 
element bearings. The program has jince undergone extensive de- 
velopment to add new capabilities and improve its execution per- 
formance. 

This user's manual describes the use of a version of the 
SHABERTH computer code developed under NASA- Lewis Research Center 
Contract NAS3-22690 with Mr. H. Coe as Technical Monitor. The 
revisions made to the program include: 

1. Modification of the cage module to calculate 
cage pocket and cage land forces i” ball, cy- 
lindrical and tapered roller bearings, using 
the models originally developed for NASA Com- 
puter Program CYBEAN [2] . An option has been 
provided to allow the specification of single 
or multiple degrees of freedom cage simulation 
by input data. 

2 . Addition of an option which permits the program 
to analyze a single ball or roller bearing 
without the specification cf shaf- geometry. 

3. Combination of two versions of code within the 
program; The SKF version and the NASA version. 

The differences between the two versions re- 
side in the calculation of the elastohydro- 
dynamic (EHD) film thickness and traction 
forces which develop }'>etween rolling element- 
raceway and rolling element-cage concentrated 
contacts. The original film thickness models 
(Archard-Cowking [3] and Dowson-Higginson [4]) 
and the Tallian traction model [5] are used 

in the SKF version, while the NASA version 
uses the Loewenthal model [6] to calculate 
film thickness and the Allen model [7] to de- 
termine traction forces. 

Additionally, a new subroutine, FLMFAC (replacing LRHS) , is 
used to determine the lubricant life factor as a function o^ A 
(film thickness/surface roughness) . The values of the lube 
life f' 7tor produced by FLMl’AC adhere closely to the curve recom- 
mended uy the ASME [8, 9). ' 
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A new section has been added. Appendix E, detailing the cal 
culation of cage pocket and cage land forces in ball and roller 
bearings. Another new section. Appendix F, describes the dif- 
ferences between the SKF and NASA methods of calculating film 
thickness and traction forces, and explains the differences in 
executing each version of the code. 
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I. INTRODUCTION 

SHABERTH performs a thermo-mechanical simulation of a load 
support system consisting of a flexible shaft supported by up to 
five (5) rolling element bearings. The shaft can be hollow or 
solid and of arbitrary geometry. Any combination of ball, 
cylindrical, or tapered roller bearings can be used to support 
the ''haft. The cylindrical roller bearing analysis permits 
thrust load to be carried by inner and outer ring flanges. The 
applied loading can consist of point or distributed moments, 
point or distributed forces, and shaft misalignments. 

Concentrated contact EHD traction models (5,10] are included 
in the program. Hydrodynamic rolling and shear forces in the 
inlet tone of the lubricated contacts are accounted for (10] . 

The effects of surface roughness (5 ], heating of the lubricant 
in the contact inlet (11], and lubricant starvation (12] are 
considered. Bearing operating clearance is determined as a 
function of shaft and housing fits, component temperatures, and 
rotational speed [13, 14] . A cage model simulates contact between 
the cage and rolling element as well as the cage and the piloting 
land. 

A lumped mass thermal model allows calculation of steady 
state or time transient system termperatures considering free 
and forced convection, conduction, radiation, and mass transport 
heat transfer [15,16,171 .A maximum of one hundred (100) temperature 
nodes can be used to describe the thermal system. 

The SHABERTH program consists of the following major 
subprograms : 

1. Bearing .Analysis - These programs are largely 
based upon the methods of Harris [18,19]. 

2. Th ree Dimensional Shaft Deflecti'.n Analysis - 
developed by Norlander and Friedrichson. 

3. Bearing Dimensional Change Analysis - based on the 
methods of Timoshenlco(13f and adapted to the shaft- 
bearing-housing system by Crecelius (14]. 

4 . Generalized Steady State and Trans i ent Temperature 
^' aprinr~and Heat Dissipation Analyses - baWd on 

z c * hods of Harris [15] Fernlund (16] and 

reason [17]. 
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2. FR0BLET1 FOR?IULATION AND SOLUTION 


The purpose of the program is to provide a tool with which 
the shaft bearing system performance characteristics can be 
determined as functions of system temperatures. These system 
temperatures may be a function of steady state operation or a 
function of time variant conditions brought on by a change 
in the system steady state condition. Such a change would be 
the termination of lubricant supply to the bearings and other 
lubricated mechanical elements. 

The program is structured with four nested, calculation 
schemes as follows: 

1. Thermal , steady state or transient temperature calcu- 
lations which predict system temperatures at a given 
operating state. 

2. Bearing dimensional eoui librium which uses the bearing 
temperatures predicted by the temperature mapping 
subprograms and the rolling element raceway load distri- 
bution, predicted by the bearing subprograms, to calcu- 
late bearing diametral clearance at a given operating 
state. 

3. Shaft-bearing system load equilibrium whiCi. calculates 
oearing inner ring positions relat'ive rc- the respective 
outer rings such that the external loading applied to the 
sha-^*: is equilibrated by the rolling element loads 
whith develop at each bearing inner ring at a given 
state. 

<1 . Bearing rolling element and cage load eouilibrium 
which calculates the rolling element ana cage equili- 
brium positions and rotational speeds based upon the 
relative inner-outer ring positions, inertia effects 
and friction conditions, which if lubricated, are 
temperature dependent. 

The above program structure allows complete mathematical 
simulation of the real physical system. The program 
has beer, coded to allow various levels of program execution 
which prove useful and economical in bearing design studies. 

These levels of execution are explained fully in Sections 
3, 4, and 5. 
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The structure of the program and the nesting of the solu- 
tion loops noted above can be seen clearly in the Program Flow 
Chart which is discussed in Appendix A. 

The sections below present the systems of field equations 
which are solved in each of the nested calculation schemes. 

A more detailed discussion is contained in [1,10 and 20]. 

2.1 Temperature Calculations 

Subsequent to each calculation of bearing generated heat 
rates, either the steady state or transient temperature mapping 
solution scheme may be executed. This sec of sequential calcu- 
lations is terminated as follows: 

1. For the steady state case; when each system temperature 
is within EPI *Centrigrade of its previously pre- 
dicted value, (EPI is specified by the user). If it is 
zero or left blank, a default value of 1° Centigrade is 
used. This criterion implies that the steady state 
equilibrium condition has been reached. 

2. The transient calculation terminates when the user 
specified time up is reached or when one of the system 
temperatures exceeds 600®C. 

2.1.j Steady State Temperature Map 

The mechanical structure to be analyzed is thought of as 
divided into a number of elements or nodes, each represented by 
a temperature. The net heat flow to node i from the surrounding 
nodes j, plus the heat generated at node i, must numerically 
equal zero* This is true fox each node i, i going from 1 co n, 
n being the number of unknown temperatures. 

After each calculation of bearing generated heat, which 
results from a solution of the shaft-bearing system portion 
of the program, a set of system temperatures is determined 
which sat.'sfy the system of equations: 

q^ ■ q^^ q • ■ 0 for all temperature nodes i (2.1) 
0 

where q . is the heat flow from all neighboring nodes to 
node i 

q i is the heat generated at node i. These values 
® may be input or calculated by the shaft bearing 
program as bearing frictional heat 
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This scheme is solved with a modified Newton -Raphson 
me hod which successfully terminates when either of two 
c:?iditions are met: 

< EP2 for all nodes i (2.2) 


where: At represents the Newton-Raphson correction to the 
temperature t at a given iteration such that, 
tN-t-i ■ tN + and N ♦ 1, and N, refer to the 
next and current iteration respectively. 

EP2 is a user specified constant. If EP2 is left blank 
or set to zero (0) a default value of 0.001 is used. 

A second convergence criterion dependent upon EP2 is also 
ustid. In the system of equations, q^j^ ♦ qgi • 0 for all nodes 
i, absolute convergence would be obtained if the right hand 
side (EQ) in fact reduced to zero (0). Usually a small residue 
remains at each node, such that (qoi q^^) ■ (EQ)i. 

The second convergence criterion is satisfied if: 



where n ■ number of equations in thermal solution 
“ number of unknown temperatures 

2 . 1,2 Transient Temperatures 

In the transient case the net heat q^^ transferred to a 
node i heats the element. It is thus necessary for heat balance 
at rode i tl ac the following equations are satisfied. 

dt. 

pi CPi V- ^ (2.4) 
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where p » density 

Cp * specific heat 

V* volume of the element 
t« temperature 
T» time 

The temperatures, t^^, at the time of initiation T • are 
assumed to be known, that is 

ti(Ts) • toi i « 1, 2, n (2.5) 


The problem of calculating the transient temperature 
distribution in a bearing arrangement thus becomes a problem 
of solving a system of non-linear differential equations of 
the first order with certain initial values given. The 
equations are non-linear since they contain terms of radiation 
and free convection, which are non-linear with temperature as 
will be shown later. The simplest and most economical way of 
solving these equations is to calculate the rate of temperature 
increase at the time T » Tj^ from equation 2.4 and then calcu^'^te 
the temperatures at time Tj^ ♦ AT from 

^ iT . tk * (2.6) 

^ pTpi: 


If the time step AT used as program input is chosen 
too large, the temperatures will oscillate, and if it is chosen 
too small the calculation will be costly. It is therefore 
desirable to choose the largest possible time step that does 
not give an oscillating solution. The program optionally cal- 
culates such a time step. The step is obtained from the 
condition, ( 16) 

^^i.k-*-! >0 i « 1, 2 n (2.7) 

°"i,k 

If this derivative were negative, the implication would 
be that the local temperature at node i has a negative effect 
on its future value. This would be tantamount to asserting 
that the hotter a region is now, the colder it will be after 
an equal time interval. An oscillating solution would result. 
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Differentiating equation (2.6) for node i, one obtains 


‘^^i.k+l * 1 + ^'^i 

at;' -- " ~ c ~r ; 


i,k 


01 pi 1 


‘^'^i k 
- — 


X^2^...n 


( 2 . 8 ) 


The derivative dq^ k calculated numerically 


dqi,k ^ qj(t^ * At^) - qi(tj) ^ 2 . 9 ' 

^^,k Ati 

For each node the value of LTi giving a value of zero to 
the right hand side of Eqn. (2.8) is calculated. The smallest 
non- zero value of AT obtained in this manner is chosen as the 
time step. 


2.1.3. Calculatior. of Heat Transfer Rate 


The transfer of heat within a medium or between two media 
can occur by conduction, convection, ratiation and fluid flow. 

All these types of heat transfer occur in a bearing appli- 
cation as the following examples show. 

1.. Heat is transferred by conduction between inner ring 
and shaft and between outer ring and housing. 

2. Heat is transferred by convection between the surface 
of the housing and the surrounding air. 

3. Heat is transferred by radiation between the shaft 
and the housing. 

4. When the bearing is lubricated and cooled by circulating 
oil, heat is transferred by fluid flow. 

Therefore, in calculating the net heat flow to a node all the 
above mentioned modes of heat transfer will be considered. 

2. 1.3.1 Generated Heat 

There may be a heat source at rode i giving rise to a heat 
flow to be added to the heat flowing from the neighboring nodes. 
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In the case that the heat source is a bearing, it may 
either be considered to produce known amounts of power, in which 
case constant numbers are entered as input to the program, or 
the shaft-bearing program may be used to calculate the bearing 
generated heat as a function of bearing temperatures. 

2 . 1 . 3 . 2 C onduction 

The heat flow q • . which is transferred by conduction 
from node i to node j , *is proportional to the difference in 
temperature (t^ - t-j) and the cross sectional area A and 
is inversely proportional to the distance i between the two 
points, thus 


‘Cl, 3 


XA 

I” 



where X* the thermal conductivity of the medium. 


( 2 . 10 ) 


2. 1.3. 3 Free Convection 

Between a solid medium such as a metallic body and a 
liquid or gas, heat transfer is by free or forced convection. 
Heat transfer by free convection is caused by the setting in 
motion of the liquid or gas as a result of a change in density 
arising from a temperature differential in the medium. With 
free convection between a solid medium and air, the heat 
energy q • • transferred between nodes i and j can be calculated 
from the^eduation, (2.11) 

q^ij - a^.A (t.-tj.) SIGNCt.-t.) (2.11) 

where * the film coefficient of heat transfer by free 
convection 

A » the surface area of contact between the media 
d ■ is an exponent, usually ■ 1.25, but any v~lue 
can be specified as input to the program 

SIGN - 
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The last factor is included to give the expression ^ 

''* » J 

a correct sign. 

The value of °^v can be calculated for various cases, see 
Jacob and Hawkins, [21] 

2. 1.3.4 Forced Convection 

Heat transfer by forced convection takes place when liquid 
or gas moves around a solid body, for example, when the liquid 
is forced to flow by means of a pump or when the solid body is 
moved through the liquid or gas. The heat flow q . . transferred 
by forced convection can be obtained from the folioAing equation. 


‘wi , j 


w 


ACt. 




( 2 . 12 ) 


where a^. is the film coefficient of heat transfer during 
forced convection. This value is dependent on 
the actual shape, the surface condition of the 
body, the difference in speed, as well as the 
properties of the liquid or gas. 

In most cases, it is possible to calculate the coefficient 
of forced convection from a general relationship of the form. 


N 


u 


aR.^P 

w 


C 


r 


(2.13) 


where a, b, and c are constants obtained from handbooks 

such as [ 22 ]. R and P are dimensionless numbers 
defined by 


Ny - Nusselt number ■ L/X 

L ■ characteristic length 
X ■ conductivity of the fluid 
R^ ■ Reynold's number « ULp/n 

U ■ characteristic speed 
p ■ density of the fluid 
n • dynamic viscosity of the fluid 
Py • Prandtl number - nCp/X 

Cp • specific heat 
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The user can input a constant value for the convection coefficient. 
Alternatively, he can let the program calculate the coefficient 
using one of the three option described below. For options 2 and 3, 
the coefficient is allowed to vary with system temperatures. 

Constant viscosity 

1 . Values of the parameters of equation (2.13) are given 
as input and a constant value of is calculated 

by the program. 

Temperature dependent viscosity 

2. The coefficient for turbulent flow and heating of 
petroleum oils is given by 

- kg . p(t)^10 (2.14) 

where Lq and kj^Q are given as input together with viscosity 
at two different temperatures. 

3. Values of the parameters of equation (2.13) are given 
as input. Viscosity is given at two different 
temperatures . 

2 . 1 , 3 . 5 Radiation 

If two flat parallel, similar surfaces are placed close 
together and have the same surface area A, the heat energy 
transferred by radiation between nodes i and j representing 
those bodies, will be, 

q^ij - coA (t^^275)^ - (t.-^273)^ (2.15) 

where c is the surface emissivity. The value of the 
coefficient c is an input variable and varies 
between 1 for a completely black surface and 0 
for an absolutely clean surface. In addition 
0 is Stefan- Boltcmann* s "radiation constant 
which has the value 5,76 x 10*8 watts/m^- (®K)^ 
and ti and t^ are the temperatures at points i 
and j : 

Heat transfer by radiation under other conditions can also 
be calculated, (21). The following equation, for instance 
applies between two concentric cylindrical surfaces. 


k 


U 
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„Ai[(t.*273)‘ * 

l"* (l-c) (Aj/aJ 


(2.16) 


where is the area of the inner cylindrical surface 
A is the area of the outer cylindrical surface 


2. 1.3. 6 Fluid Flow 


Between nodes established in fluids, heat is transferred 
by transport of the fluid itself and the heat it contains. 

Qoi 

♦ - - ir 

I ©i 0j i 

r—7—/ — / V / — 7 — 7 " 

FIGURE 2.1 CONVECTIVE HEAT TRANSFER 


Figure 2.1 shows nodes i and j at the midpoints of consecu- 
tive segments established in a stream of flowing fluid. 

The heat flow through the boundary between nodes 

i and j can be calculated as the sum of the heat flow q£^ through 

the middle of the element i, and half the heat flow q^. trans- 

^01 

ferred to node i by other means, such as convection. 

The heat carried by mass flow is, 

'^fi • Si S (2.17) 

where Vj ■ the volume flow rate through node i 


(1 

11 

(I 
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The heat input to node i is the sun of the heat cenerated 
at node i (if any) and the sun over all other nodes of the 
heat transferreo to node i by conduction, radiation, free and 
forced convection. 


*oi 


Og, 


• r (q • • 
1 >1 






(2.18) 


The heat flow between the nodes of Fig. 2.1 is then 


‘ui,j 


qfi 


q -/2 


(2.19) 


If the flow fron node i is dividing between nodes j & k,(Fig. 2.2) 
then the heat flow is calculated fron 




•'ij 


(qfi 


• Ooi/2) 


C2.20) 


where K. ■ the proportion of the flow at i going to 
node j, 0 <K. . <1. K. . is specified at 

input. •' " 



2. 1.3. 7 Total Heat Transferred 


The net heat flow rate to node i can be expressed as. 


*^i * ^G,i * * ^ui,j * ^vi,j *^wi,j * *^Ri,)) (2»21) 
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The summation should include all nodes j, both with 
unknown temperatures as well as boundary nodes, at which the 
temperature is knowi\ so long as they have a direct heat exchange 
with node i. 

This expression is a non-linear function of temperatures 
because of the terms qy^ and q„. Therefore the equations to be 
solved for a steady state solution are non-linear. The sub- 
program SOLITCX for solving non-linear simultaneous equations is 
used for this purpose. 

2.1.4 Conduction Through a Bearing 

As described in Section 2. 1,3. 2 the conduction between 
two nodes is governed by the thermal conductivity parameter \ 
of the medium through which conduction takes place. The value 
of X is specified at input. 

An exception is when one of the nodes represents a 
bearing ring and the other a set of rolling elements. In this 
case the conduction is separately calculated using the princi- 
ples described below. Note that separate calculations are 
performed for the rolling element raceway contacts and the roll- 
ing element-flange contacts. The methods for both calculations 
are identical and are performed within the program. 

2. 1.4.1. Thermal Resistance 


It is assumed that the rolling speeds of the rolling 
elements are so high that the bulk temperature of the rolling 
elements is the same at both the inner and outer races, except 
in a volume close to the surface. The resistance to heat flow 
can then be calculated as the sum of the resistance across 
the surface and the resistance of the material close to the 
surface . 


The resistance is defined implicitly by 

dt«fl»q (2.22) 


where 

Lt is temperature difference 
q is heat flow 

The resistance due to conduction through the EHD film is 
calculated as 


(2.23) 


- h 

TT 
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where h is taken to be the calculated plateau film thickness 
A is the Hertzian contact area at the specific rolling 
elementTing ontact under consideration. 

X is the conductivity of the oil. 

The geometry is shown in Figure 2.3(a). 


So far, a constant temperature difference between the 
surfaces has been assumed. But during the time period of 
contact, the difference will decrease because of the finite 
thermal difrusivity of the material near the s'uface. Fig. 2.3 
(b). 


To points at a distance from the surface, this phenomenon 
will have the same effect as an additional resistance ^2 acting 
in series with fli. 

This resistance was estimated in [23] as. 


^ 2 “ ^ 



h 


( 2.2 


A > 


where 

X 

P 


contact length, or in the case of an 
elliptical contact area, 0.8 times the 
major axis 

heat conductivity 

thermal diffusivity « X/(p.Cp) 

density 


Cp ■ specific heat 
b ■ half the contact width 


V ■ rolling speed 


The resultant resistance is 


n 


res 




a. 


(2.25) 
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There is one such resistance at each rolling element. 
They all act in parallel. The ®Q*^ivalent resistance eqv 
is thus obtained from 



n 

E 1 
i-1 "TJres.i 


(2.26) 


2.2 Bearing Dimensional Change Analysis 

The program calculates the changes ir bearing diametral 
clearances according to the analysis described in [14], and 
expressed in generalized equation form as, 

ADCL - f {(Fits)^, ^i»p. m» ^Vm^» ® (2.27) 

and outer 

rings respectively 

i • 1,2, 3, 4, 5 for 
shaft, inner ring, 
outer ring, housing 
and rolling element 
respectively 
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where: i(iDCL is the change in bearing diametral clearance 

Fits are the cold mounted shaft and housing fits, 
tj^ are the component temperatures 

refers to the ring rotational speeds 
Qy refers to the radial component of the minimum 
rolling element-race normal force 


A bearing clearance change criterion is satisfied when 
the change in bearing diametral clearance remains within 
a narrow, user specified range, for two successive iterations 
as follows: 

I (riDCL)j^ - CADCL)j^.jj< EPSFIT for all bearings (2.28) 


D 

where: N denotes the most recent iteration and 

N-1 denotes the previous iteration, 

D denotes the ball or roller diameter and 
EPSFIT is a user specified value, (default value • .OOOID) 

It should be noted that although ring rotational speeds, 
and initial, i.e. cold, shaft and housing fits are considered 
in the clearance change analysis, these two factors are fixed 
at input and remain constant through the entire solution. 

Although component temperatures may change as a consequence 
oi the thermal solution, temperatures remain constant through 
a complete set of clearance change iterations. As a result, only 
the change in bearing load distribution affects the change in 
bearing clearance within a set of clearance change iterations. 

“ • ^ Bearing Inner Ring Equilibrium 

The bearing inner ring equilibrium solution is obtained 
by solving the system: 

- (FM^)i ■ 0 for all bearings, i (2.29) 

where: FM^, denotes a vector of bearing loads and 

moments resulting from rolling element/race 
forces and moments. 
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^bi • 


bxi 


byi 

c 

bzi 


M 


byi 


M, 


bzi 


Forces 


Moments 


( 2 , 30 ) 


If the bearing solution considers friction, FMv is com- 
prised of the ball race friction forces as well as the 
normal forces. 


If the beating solution is, at the user’s option, 
frictionless, FM. is comprised only of rolling element/race 
normal contact forces. 


FM i denotes a similar vector of loads, exerted on the 
inner ring by the shaft. 


™si 


sxi 


syi 


szi 


M. 


syi 


M 


szi 


Forces 


Moments 


(2.31) 


The variables in this system of equations are thc^bearing 
inner ring deflections d. and the shaft displacements dg at 
all bearing locations. The bearing loads may be expressed as 
a function of the inner ring deflections. 

FMj^ - FMjj (d^) (2.32) 
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The deflection a bearing is described by two radial 

deflections and 5 2 angular deflections 6 y and 6j and an 
axial deflection axial deflection is assumed to oe 

the same for all bearings and the shaft. 


The solution scheme is ended when 


(''’ij < 

■ 


If for some i or 
( 2.33) . 


EPSl f frictionless) 
EPS2 (friction) 


(2.33 ) 


i ■ 1,... (Number of bearings) 

j ■ 1,5 - for the 3 linear and two 
angular deflections at 
each bearing 

j, (d) • 0, Eq (2.34) is used in place of 


^ ij ^EPSl (frictionless) 
(O.OCl X NBR^ ESP2 (friction) 


NBRG denotes the number of bearings in the system. 


EPSl or EPS2 is used depending on whether the bearing solu- 
tions are fricitonless or include friction, respectively. If 
the bearing deflections are extremily small, computer-generated 
numerical inaccuracies may prevent convergence according to the 
above criteria although a perfectly good solution has been ob- 
tained. To overcome this problem, the iteration is terminated 
if all angular deflections are less than 2 x 10-^ radians 
and all linear deflections are less than 5 x 10-® inches. Any 
one of the above criteria imply that inner ring equilibrium is 
satisfied. 

- • ^ B earing Quas j-Dynamic Solution 

The bearing quasi-dynamic solution is obtained through a 
two step process: 

1) Elastic Solution - considering rolling element centri- 
fugal force, plus the gyroscopic moment for a tapered 
roller, 

2) Elastic and Quasi -dYTiamic Solution* 

"Quasi -dynamic equilibrium is used to connote that tne true 
dynamic equilibrium terms containing first derivatives of the 
ball rotational speed vectors and the second derivatives of 
rolling element position vectors with respect to time are 
replaced by numerical expressions which are position rather 
than time dependent. 
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The equations which define rolling element quasi-dynamic 
force equilibrium takes the form 


cn 

n 


(Q. 


m 


where :Q 




m » 1-3 refers to the outer 
inner and cage rolling element 
contacts respectively 
m ■ 1-4 for tapered rollers with 
one roller end- flange contact. 

+ F_ l+F^O ni ■ 1-7 for cylindrical rollers 

\'here four roller and flange (2.35) 

contacts are possible. 


■.} 


m 


m 


is the vector normal load per unit length 
of the contact. See Ref. [ij. 

is the vector of friction force per unit length 
of the contact. See Ref. [10]. 


F is the vector of inertia and drag forces. See Ref,[l) 
s is a coordinate along the contact 

perpendicular to the direction of rolling (usually 
the major axis) 


^ is half the contact length. See Ref. II]. 

Fjjj is the vector sum of the hydrodynamic forces 

acting on the rolling element at the m-th contact. 
For ball-raceway contact see Ref. 110]. For the 
roller- raceway contact, see Ref. [21] . 


Rolling element moment equil’^rium is defined by: 


r 

m 




♦ r_ X F_ ♦ M, 

m m I 


(2.36) 


Q_ f_ F a and s are defined above, Mr is a vector 

* * * ' of inertia moments. For the definition of Mj 

... refer to Ref. [ l] . 

r is a vector from the rolling element 
center to the point of contact. 


In the frictionless elastic solution F and f ■ 0, 

m m 

Additionally, the only rolling element inertia term considered 
in the frictionless solution is centrifugal force, plus the gyro- 
scopic moment for tapered rollers. As a consequence only the 
axial and radial force equilibrium equations are solved for each 
ball. For each roller the radial and axial force equilibrium 
and the tilting moment about the t axis of Fig. 2.4 is solved. 
A dummy equation for axial force equilibrium is included in the 
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FIGURE 2 .^ 

BEARING INERTIAL (XVZ) AND ROLLING ELEMENT (xyz) , 

COORDINATE SYSTEMS 
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solution matrix which keeps the roller centered with respect to 
the outer race if the cylindrical roller bearing takes no thrust 
load. 

The friction solution determines ball quasi-dynamic equilib- 
rium for six degrees of freedom. The rolling-element variables 
in this solution are xl, yl, wx, wy, toz, and wo. 

where x.is the rolling element axial position relative 
to the outer race groove curvature center, 
y.is the rolling element radial position relative 
to the outer race groove curvature center. 
u)X, (ijy, toz are orthogonal rolling element rotational 
speeds relative to the cage speed, about the x, y, and 
z axes and Wo is the rolling element orbital speed. 

The variables x. and y^ are the ball unknowns in the frictionless 
solution. The Variables in the roller frictionless solution are 
*1' ^1' ~ arctan (coy/wx) 

Details of the cage analysis are contained in Appendix E. 

Either one or three cage equilibrium equations are considered, de- 
pending upon the number of degrees of freedom given to the cage. 

The cage equations and cage-rolling element interactions are 
not considered when the friction forces are omitted from the rol- 
ling element equilibrium equations. 

2.4.1. Cage Degrees of Freedom 

The program has been modified to allow the user to specify the 
number of degrees of freedom (DOF) of the cage as either 1 or 3 . 
The single degree of freedom corresponds to a smaller angular rota- 
tion about the bearing axis, measured with respect to rolling 
element 1. The angular displacement is converted to a linear di- 
mension by multiplying it by the bearing pitch diameter and is 
noted in Fig. 2.4 as z . When a single degree of freedom is in- 
put, the sum of momentS acting on the cage about the bearing X 
axis is required to be zero. This moment equation considers the 
cage-rolling element normal and friction forces as well as the 
torque generated at the cage-ring surface. 

If the user assigns three degrees of freedom to the cage, it 
is permitted to move to an eccentric position with respect to the 
land on viuch it is piloted. The additional degrees of freedom are 
the cage center of mass radial displacement, e, and the angular 
displacement 0 of the center of mass, with respect to the bearing 
Y axis. (See figures 2.5 and 2.6.) The radial friction forces 
as well as the pressure build-up between the cage and its piloting 
surface are considered in the equilibrium equations. The effect of 
the cage mass is neglected. 
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Occasionally » the program will determine that an input value 
of three degrees of freedom is impractical, in which case it will 
override the input and allow the cage to have only one DOF. This 
occurs if the cage will tend to rotate concentrically with respect 
to the ring on which it is riding. Such a condition is determined 
as a function of the rolling element orbital speed variation and 
prevails with most roller bearings and with ball bearings subjected 
only to axial loading. In both cases, orbital speed variation is 
often inconsequential. Also, when the cage rides on the rolling 
elements, it is allowed only one degree of freedom. 

The ball bearing friction solution is thus obtained by solving 
6Z+(1 or 3) equations where Z is the number of rolling elements. 

The ball bearing frictionless solution is obtained by solving 1, 
(Z/2) (Z/2-i-l) or Z sets of 2 equations, depending upon the number 
of rolling elements in the bearing and the degree of load symmetry 
which prevails. The various symmetry conditions are explained be- 
low. 


The roller bearing friction solution contains 4Z+(1 or 3) 
equations and the frictionless solution contains Z/2, Z/2+1 or Z 
sets of three equations again depending upon the number or rolling 
elements and whether or not load symmetry exists. 

The various load symmetry conditions are as follows. Axial 
symmetry is utilized if the load is axial only, then only one set 
of two or three equations is solved for the frictionless case and 
six rolling element and one cage equilibrium equations are solved 
when friction is included. All rolling elements are assumed to 
behave identically. 

Radial load symmetry is utilized if the non-axial shaft load- 
ing is comprised of only radial components parallel to the Y axis 
and moment components parallel to the Z axis and the position of 
the first rolling element is utilized. When this symmetry exists, 
only half the rolling elements need be considered if the number 
of rolling elements is even and one half plus one need be considered 
if the number is odd. Because of inertia terms, radial load sym- 
metry can only be utilized in the frictionless calculations. 

If load symmetry is not present, then Z sets of two (ball 
bearing) or Z sets of three (roller bearing) equations must be 
solved to obtain the frictionless solution. 

As with the steady state temperature mapping scheme, the 
Newton-Raphson scheme in subprogram S0LV13 is used to solve the 
sets of equations for each bearing. The iteration scheme term- 
inates when either: 
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IKii 


(K-1)^ 


EPSl frictlonleas 


or 


^ 2 
I EQ‘ 

Ii=l 


L n j 


(2.37) 


i=l...n 

.EPS2 friction 
Where K and K-1 refer to iteration numbers 


H 

< 


100 * 


C 


EPSl frictionless 


EPS2 friction 


(2.38) 


Experience has shown that the ser;ond criterion is usually 
responsible for terminating the solution. However, when rolling 
element loads are extremely large, on the order of 10^ Newtons, 
it becomes difficult to reduce the equation residues to less 
than 10 Newtons. In those instances, the first criterion usually 
terminates the iteration scheme. 
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3.1 Types of Input Data 

A complete set of input data comprises data of four distinct 
categories. Within these categories , cards which convey specific 
kinds of information are referred to as card types. Depending 
on the complexity of the problem./ the input data set may contain 
none, one or several cards of a given type. The categories are 
listed below. 

I . Title Cards 

A title card plus a second card which provides the 
program control information for the shaft-bearing 
solution. 

II . Bearing Data Cards 

A set of up to sixteeen (16) card types. Each set 
describes one bearing in the assembly. All bearings 
must be so described. The card sets must be input 
sequentially in order of increasing distance from a 
selected end of the shaft. 

Ill . Thermal Data Cards 

A set of up to nine (9) card types to describe the 
thermal model of the assembly. 

IV. Loading Data Card 

One card describing the loading on a single bearing. 

This card type is used only when Shaft Data Cards are 
omitted . 

V. Shaft Data Cards 

A set of three (3) card types to describe the shaft 
geometry, bearing locations on the shaft and shaft load- 
ing. Used only when Loading Data Card is omitted . 

If the program is being used to predict the performance of a 
bearing assembly, cards from sets I, II, III, and V must be in- 
cluded in the runstream. If the program is being used to thermally 
model a mechanical system wherein no bearing calculations need be 
performed, the cards from sets II, IV and V are omitted. 

The review of required input information which follows is 
broken into the five sets of data categories given above, with 
special emphasis on program control data. 

The input data instructions are given in Appendix C, and are 
for the most part, self-explanatory. They are laid out in the 
format of an eighty column data card . A description of the var- 
iables is given in the input instruction forms. 
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The units used for input data are as follows: 

Linear Dimensions - (mm) 

Angles - (degrees) 

Surface Roughness (microns) 

Bearing Angular Mounting Errors - (radians) 
Rotational Speeds - (RPM) 

Force - (Newtons) (N) 

Moments - (N-mm) 

Pressure, Elastic Modulus - (N/ram^) 

Density - (gm/cm^) 

Kinematic Viscosity - (cs) 

Temperature - (degrees centigrade) (®C) 
Coefficient of Thermal Expansion - (OC*-) 
Thermal Conductivity - (Watts/m/®C) 

3.2 Data Set I - Title Cards 


3.2.1 Title Card 1 


This card should contain the computer run title and any 
information which might prove u.-eful for future identification. 
The full eighty (80) columns are available for this purpose. 

The title will appear at the top of each page of Program output, 

3.2.2 Title Card 2 


This card provides the control infonnation for the shaft 
bearing solution. 

Item 1: Shaft Speed in rpm, GOV (1). All bearings have 

the same sEaTt and inner ring speed. 

Item 2 ; Number of Bearings on the shaft (NBRG), a minimum 
of zero is permitted if no bearing solution is being sought. A 
maximum of five is permitted. Note that a bearing is defined as 
a single row of rolling elements. Thus a double row bearing is 
treated as two separate, single row bearings. 

Item 3: Print Flag (NPRINT) , NPRINT equal to zero is normal 

and will result in no intermediate or debug output. With a value 
of one, a low level intermediate print is obtained at the end cf 
each shaft bearing iteration. The values of the i-nner ring dis- 
placements (DEL), equation residues, bearing inner ring residual 
loads, and bearing partial derivatives are printed for each iteration. 
This level is recommended in cases where proper convergence of the 
bearing solution does not occur. 

At the end of each bearing iteration, wherein the roiling 
element and cage equilibrium equations are solved, an error para- 
meter is printed which has the value: 
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Error Parameter • 

is the change in the variable X specified at 
iteration N. 

%-l value of the variable specified at 

the previous iteration. 

The Error Parameter is calculated for each of the bearing 
variables, but only the largest one is printed. 

Additionally, at the end of each Clearance Change itera- 
tion, the clearance change error parameter is printed. This 
error is defined: 


Error Parameter 


Rolling Element Diameter 


where DCLjj and DCLjyj.j denote the clearance changes calculated 
at the current and previous iterations respectively. 

If NPRINT is set at 2 all of the above information is 
printed. Additionally the variable values and residue values 
are printed for each iteration of the rolling element and cage 
equilibrium solution. This level is not generally reconanended be- 
cause of the large volume of output produced. 


Item 4 : ITFIT controls the number of iterations allowed to 
satisfy the bearing clearance change iteration scheme. If 
ITFIT is set to zero (0), or left blank, the clearance change 
portion of the program is not executed. If a positive integer 
is input, the clearance change scheme is utilized with a maximum 
iteration limit of five (5). If a negative integer is input, 
the scheme is used with a maximum iteration limit equal to 
the absolute value of the negative integer. 

Item 5 : ITMAIN limits the number of iterations attempted 
during 'tTie solution of the shaft and bearing inner ring equili- 
brium problems, i.e. establishing the equilibrium of bearing 
reactions and applied shaft loads. If ITMAIN is left blank, 
set to zero, or to a positive integer, then (15) iterations 
are permitted. If ITMAIN is set to a negative integer the 
number of iterations is limited to the absolute value of that 
integer. 


30 



ATS IDO 40 


Item 6 ; GOV (2) or EPSFIT is the convergence criterion for 
the diametral clearance change portion of the analysis. As men- 
tioned under item 3 above, this error parameter is defined by Eq. 
2.28. 


The iteration scheme is terminated when the error parameter 
is less than the input value of EPSFIT. If EPSFIT is left blank 
or is set to zero (0), the program default value of 0.0001 times 
the rolling element diameter is used. 

Items 7 & 8 ; Main loop accuracy for frictionless elastic 
(EPSll and friction solution (EPS2) . These accuracy values con- 
trol the accuracy of the shaft bearing deflection solution as 
well as the quasi-dynamic solution of the component dynamics. If 
EPSl and EPS2 are left blank or set to zero (0) , default values 
of 0.001 and 0.0001 respectively are used. 

Item 9 ; JUSTBR, column 78, is a flag indicating whether or 
not a single ball or roller bearing is to be analyzed. If a 
value of 1 is input, data are needed for only one bearing, and 
loads are input on Loading Data card Ll following the Thermal 
Data cards. Also, Shaft Data cards are omitted. If JUSTBR is 
0 or blank, up to 5 bearings are analyzed with Shaft Data, and 
the Loading Data card is omitted. 

Item 10 ; IMT, if set to 1, the material properties for both 
bearing rings and the rolling elements are to be input on card 
types B 11 through B 14. If IMT is zero or blank, the rings and 
rolling elements are assumed to be 52100 bearing steel. Card 
types B 11 through B 14 are required if the change in bearing 
diametral clearance is to be calculated or if a system component 
has properties different from steel. 

Item 11 ; NPASS controls the level of the bearing solution: 

0 Elastic Contact Forces are calculated. No lubrication 
or friction effects are considered. 

1 Elastic Contact Forces are calculated. Lubrication 
and friction effects are considered using raceway con- 
trol (ball bearing) or epicyclic (roller bearing) 
assumptions to estimate rolling element and cage speeds. 
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2 Inner Equilibrium is satisfied considering only 
the Elastic Contact Forces. Using the inner ring 
positions thus obtained, rolling element and cage 
equilibrium are determined considering friction. 

3 Complete Solution. The inner ring, rolling element 
and cage equilibrium is determined considering all 
elastic and friction forces. 

3.3 Data Set II - Bearing Data 

Most of the input instructicas are self-explanatory. Where 
certain items are deemed to require more explanation than given 
in the input data format instructions they are treated on an in- 
dividual basis by card type and item number. 

Most of the bearing input data is read into a two dimensional 

array named "BD," which has the dimensions (1830, 5). For each 

of the five bearings permitted on a shaft, a total of 1830 pieces 

of data may be stored. Denoting BD(I,J), I represents a specific 

piece of bearing data, J represents the bearing number. The 
bearing input data of Data Set II occupies the first 106 locations 
of the 1830 allotted. On the input data format sheets the des- 
ignation BD(I) where 1=1... 106, denotes the location within the BD 
array where each piece of input data is stored. 
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3.3.1 Card Tv-pe 1 - Bearing Type and Material Designations 

Item 1 : Bearing type, columns 1-10 must be specified, 

left justified, i.e., "B”, ”C" or *'T*' in column 1. This format 
must be followed since the Program recognition of bearing type, 
(ball, cylindrical or tapered roller bearing), is derived from 
reading the ”B", "C”, or "T" in the first column of this card. 

Items 2 S 3 : Columns 11-30 and 31-50, ’’Steel designations”, 
inner and outer rings respectively. The alphameric- literal 
description of the steel types such as "M-50” or ”AISI 52100” 
is input. 

Items 4 S 5 : Columns 51-60 and 61-70, the numbers input 
for items 4 and 5 are used to account for improved materials 
and multiply the raceway fatigue lives as determined by 
Lundberg-Palragren methods. Typical life factor values for modern 
steels are in the neighborhood of 2.0 to 3.0. If the ASM£ 
Publication Life Adjustment Factors for Ball and Roller Bearings , 
is reference?n>y the user, tHe ^iaterial J^actor D and the 
Material Process Factor E should be used multiplicatively 
as inputs for items 4 and 3. The program computes a lubricant 
life factor based on the value of A (EHD plateau film thickness/ 
composite RMS surface roughness). The calculated lube life 
factor ranges from 0.21 for A < 0.6 to 3.0 for A > 10.0. 


Item 6 ; Columns 71-78, "Orientation angle of the first 
rolling element". (4>1) (degrees). Refer to Fig. (2.4). The 
quasi-dynamic rolling element bearing problem has an infinite 
number of solutions which fall within a narrow envelope having 
a periodic shape. The solution obtained is a function of 
the rolling element positions relative to the bearing system 
coordinate axes. 41 ■ 0 places a rolling element on the Y 

axis and is the choice customarily made. 4 ^ can be desig- 
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nated as any value 0<4ii< 360/Z where Z is the nuinber of rolling 
elements. For each JitYerent value assigned to a different 
ent, although similar, bearing solution will be obtained. To 
take advantage of bearing symmetry and the computer time saving s 
which result, must be specified as zero or lett blank . 

Item: 7 Column 80, a signal, termed the crown drop fls.g, 

which specifies for a cylindrical or tapered roller bearing, 
whether the roller-race crown drops will be calculated, or read 
directly. If item 7 is blank or zero, the crown drops are cal- 
culated based on the roller-race crown radius, and effective flat 
length input information. If the crown drop flag is other than 
zero or blank the non-uniform separation of the roller and race- 
way must be specified at the center of each slice into which the 
roller-raceway effective contact length is divided. The slice 
widths are identical. The number of slices is input as item 7 
card type B4. The non-uniform roller-raceway separation is in- 
put on card types B5 and B6, 

3.3.2 Card Type B2 - Bearing Geometry and Outer Ring Speed 
3. 3. 2.1. Ball Bearing Geometry 

Items 1, 2 and 6 are self explanatory. Item 5 pertains 
only to a tapered roller bearing, as discussed later. Items 3 
and 4 require explanation however. 

Through the proper specification of the diametral clearance 
and contact angle, the Program can properly handle deep groove, 
split inner, and angular contact ball bearings. 

The deep groove ball bearing requires the specification 
of the contact angle corresponding to either the operatino diame- 
tral clearance Pd or the off the shelf diametral clearance, if the 
dimensional change analysis is utilized. 

The angular contact bearing is fully described through 
specification of the contact angle which obtains under a gauge, 
axial load. However, this method of input does not accurately 
define the system if there is more than one angular contact 
supporting the shaft and at least one of those bearings has 
its grooves offset in the direction opposite to the other 
bearings and if the shaft is capable of axial and/or radial 
play. In other words, if, what are known as angular contact 
ball bearings, are mounted such that some diametral shaft 
play is permitted, an auxilliary angle as well as the diametral 
play must be specified at input. The angle input is not the 
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manufacturer’s designated contact angle, ot , but an auxilliary 
angle, the calculation for which shall be demonstrated. 

Refer to Figure 3.1. The manufacturer’s contact angle is 


i calculated as follows: 


[■ 

^ ..fZA-'Pi'l 
o<. s Cq& L ZA j 

C3.1) 

r~ 

1. 

A. a To + ri - D 

C3.2) 

f 

[ 

where: rQ and r^ are the outer and inner 

radii respectively 

racew.y groove 

f" , 
! 

D is the ball diameter 


l 

Under a gauge axial load, o is obtained at both inner and 
outer raceways for each ball. Under this condition, the outer 
and inner raceways are axially offset an amount S® . 

■ - 


(3.3) 


\ 


i 


When angular contact ball bearings are mounted with some 
diametral play, the grooves are offset an amount Soq such that 
Sao Q* The diametral play which obtains at this condition 
is Sd. This diametral play is usually known by the engineer 
or designer and is usually required to allow some forgiveness 
when thermal gradients are encountered. Assuming that the user 
has the values for a , rQ, r^, D and Sqq then: 



where: Pd and A may be calculated from Eqs. (3.1) and (3.2). 


[ 

I 

I 


If Stto is unknown, the following equation may be solved 
for CIO- 
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FIGURE 3.1 ANGULAR CONTACT BALL BEARING GEOMETRY 





Geometry Under a Gauge Axial 
Load Showing the Manufacturer’s 
Contact Angle. 


Geometry with Grooves Aligned 
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FIGURE 3.1 ANGULAR CONTACT BALL BEARING GEOMETRY 
(CONTINUED) 




(c) 

Geometry Showing the as Mounted 
Groove Alignment* the Input 
Diametral Play (Sd), and the 
Input Contact Angle (a©)* 
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a 


o 


cos“l [ 


2A-P<1 + Sd 1 

— n J 


(3.5) 


In order that the Program properly handle split inner ring 
ball bearings an auxiliary angle and diametral play must be 
input. Referring to Figure 3.2, the auxiliary angle Oq and 
diametral play Sd must be determined and input. Typically the 
values of D, r^, r^, as and Sd' (assembled bearing diametral play) 
are known. 


The unloaded half of the inner ring must be removed from con- 
sideration and the ball moved such that its center lies on the 
line connecting the origins of ri and ro and positioned such 
that the auxiliary clearance Sd/M exists at both the inner and 
outer raceways. The auxiliary angle is given by: 


- tin”l r(ri“D/2)sin Cs 

- ^an L= — *n7 T- g;i TV' » T~' ■ 


] 


ro-D72^Sd»7Z + (ri-D72)cds Os ^ (3.6) 


The input bearing diametral play, Sd , can then be calculated 
as follows: 


Sd = Sd' + (2ri-D)(l-cos Og) - 2A(l-cos «o) (3-7) 


The angle associated with each ball bearing must be 
specified with the correct sign. A positive contact angle 
allows the bearing to accept a positively directed axial load 
transmitted by the shaft. 


SKF TECHN(3L(3GY SEPVCES 
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FIGURE 3.2 SPLIT INNER RING BALL BEARING GEOMETRY 
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3.3. 2.2 Tapered Roller Geometry 

Items 2 and 6 are self explanatory. 

Item 1 ; Bearing Pitch Diameter 

The tapered roller bearing pitch diameter may be calculated 
through specification of the roller measured large end diameter » 
the roller large end corner radius » one half the cup included 
angle, the roller included angle and the roller total length. 

This calculation is performed within the program. The user there- 
fore need not specify the tapered roller bearing pitch diameter. 

Item 3 : Bearing Axial Play 

For a tapered roller bearing, the bearing axial play rather 
than diametral clearance must be specified. (Sec Fig. 3.3) Note 
that this end play pertains to only the bearing in question. For 
two identical tapered roller bearings on a shaft, one half the 
total shaft axial play should be specified for each bearing. For 
two tapered roller bearings of dissimilar site, the tota . shaft 
axial play should be apportioned according to the bearing site 
such that the sura of the axial plays specified for the two 
bearings equals the total for the shaft. 

Item 4 : Bearing Contact Angle 

For the tapered roller bearing, one half the incliided cun 
angle (o ) is input as the contact angle. (S^e Fig. 5.o) 

This ancle must be specified with the correct sign. A positive 
angle allows the bearing to accept a positively directed axial*load 
transmitted by the shaft and vice versa for a negative angleT 

Item 5 : Tapered Roller Bearing Flange Angle 

The flange angle is shown by o f in Fig. 3.3. The flange 
angle must always be positive. 


3 . 3 . 2 . 3 Cylindrical Roller Bearing Geometry 

Items 1, 2, 3 and 6 are self explanatory. Both items 4 
and 5 should be left blank. 
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FIGURE 3.3 TAPERED ROLLER BEARING GEOMETRY 
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3.3.3 Card Tvpe B3 - Rolling Element Geometr) 


3. 3. 3.1 Ball Geometry 

The geometry of a ball is fully defined by its diameter. 

3. 3. 3. 2 Tapered Roller Geometry 

Item 1 - Roller. Measured Large End Diameter 

For a tapered roller, the input diameter required is the 
largest measureable diameter. (See Fig. 3.4.) Typically 
this measurement should be taken where the large end corner 
radius becomes tangent to the roller surface profile. Within 
the program a "working*' large end diameter is calculated. 

This diameter is shown as in Fig. 3.4. All bearing geometrical 
relationships are calculatea based on Dw- 




where : 


(3.8) 


D is the measured large end diameter, reo is 
the distance from the roller end to the 
beginning of the roller effective length at 
the outer raceway surface, measured parallel 
to the roller surface. If r© is the corner 
radius at the roller large end 


- 




(3.9) 


V is one half the roller included angle. 


Item 2 through 
e xp 1 anatory^ Note 
flat length, Item 6, 
crowned. 


through 6 are shown in Fig. 3.4 and are self 

Note that the program can handle a nonzero roller 
Item 6, Most tapered rollers are, however, fully 
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3. 3. 3. 3 Cylindrical Roller Geometry 

Typically, cylindrical rollers are partially crowned as 
:;hown in Fig. 3.5. The center of the roller is flat. Toward 
the ends the roller profile is formed by a crown radius, R . 

There are usually rounded corners at the roller ends. These 
corners reduce the load carrying surface of the roller such that 
if there are no raceway undercuts the roller raceway effective 
length equals the roller ^otal length less the two corner radii 
(see Section 3. 3. 4. 2). Note that a partially crowned roller is 
specified through input of a non-zero flat length. If the flat 
length is zero, the roller is fully crowned with its profile 
defined by the roller crown radius. 

The roller end sphere radius, R , is defined in Fig. 3.5. 

For a cylindrical roller, the roller®included angle is zero. 

3.3.4 Card Type B-4 - Rolling Element-Ring Geometry 

3. 3. 4.1 Ball Bearing 

Items 1 and 2 refer to the outer and inner raceway curva- 
ture s“reipecFIveTy where curvature is defined as the cross 
groove radius divided by the ball diameter. Typical values 
range from 0.515 to 0.57. 

3. 3. 4. 2 Tapered Roller Bearing Contact Geometry 

Items 1 and 2 - Roller Raceway Effective length ^ 

The roller-raceway load bearing surface is measured parallel 
to the roller surface such that if there were no relief at the 
roller ends the effective contact length would be: 

L| = L/cosv (3.10) 

However, since the roller has corners at the large and small 
ends, the actual effective length is less than L§. The con- 
sideration of raceway undercuts at the inner raceway flanges 
may result in a raceway effective length less than the roller 
effective length in which case the shorter of the two should 
be input. 
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FIG. 3.5 CYLINDRICAL ROLLER BEARING GEOMETRY 
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Items 3 and 4 - Raceway Crown Radius 

The present analysis permits both the roller and race- 
way to be crowned. If the raceway is crowned, it must be fully 
crowned with no flat specified. If the raceways are flat the 
input crown radius may be left blank, in which case a default 
value of 1. X lO'*’ inches is used. 

Note that the unloaded roller-raceway separation along 
the roller profile, (5c) Fig. 3.6, calculated at the center of 
each roller raceway slice is comprised of the sum of the roller 
and raceway crown drops. 


Items 5 and 6 - Roller Large End Corner Relief 

These data specify the distance from the roller large 
end to the point on the roller surface where the roller effective 
length begins. For the outer raceway contact this distance 
inav be calculated using Eq . 3.9. For the inner raceway use Eq.3.9 
or the width of the inner raceway undercut at the large end. (Fig. 3.4) 

Item 7 - The number of slices into which the roller 
raceway contacts are divided. 

A maximum value of twenty (20) is permitted. A default 
value of eleven (11) is used if a blank or zero is read. 

3. 3. 4, 3 Cylindrical Roller Contact Geometry 

3. 3.4. 3.1 Card Type B4A Roller Raceway Geometry 

Items 1 through 4 and 7 have the same definitions as 
they had for the tapered roller bearing. 

Items 5 and 6 - Roller End Corner Relief is not re- 
v^uired innut for the cylindrical roller bearing since the roller 
raceway effective length is assumed to be centered along the 
roller. This was not the case for the tapered roller bearing. 

Card Type B4B Roller Flange Geometry for 
Cylindrical Roller Bearings 


3.3.4. 3. 2 
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Items 1 - 4 - The flange angle and end play definitions 
can be seen in Fig. 3.5. Note that when a ring has a single 
flange such as Inversion No. 3. Fig. 3.7, the end play is the 
distance between the roller and the flange when the roller is 
centered on the raceway. 

Item 5 • The Flange Inversion Index - See Fig. 3.7 
The number which corresponds to the particular Flange inversion 
being examined must be input. Note that since the inversions 
greater than eight (8) cannot carry axial loading, the bearing 
carries load only on the raceways, and thus the Program resets 
the inversion index to one if the input value is greater than 
eight. The inversion index must be input as a real number, (with 
a decimal point j . 

3.3.5 Roller-Raceway Non-Uniform Profile Definition 
3. 3. 5.1 Card Type B5 - Outer Raceway Roller Contact 

These cards are used to input the separation between the cute 
raceway § roller at the center of each slice along the roller profile 
with the high points of the roller and race in contact, i.e., with 
all clearance between roller and raceway removed. These cards 
must be omitted if item 7 of the Bearing Data Title card is zero 
or blank. These data are used only when the roller-raceway 
profile geometry cannot be defined by card types B5 and b4a . 

3.3.6 Card Type B6 - Inner Raceway Roller Contact 

Same as Card Type B5 for the inner raceway-roller contact. 


3.3.7 Ring-Rolling Element Surface Data 

5. 5. 7.1 Card Type B7A - Raceway - Rolling Element Surface Data 

Items 1 through 6 define the statistical surface micro- 
geometry parameters of the rollers and raceways. Items 1 
through 3 require the input of center line average CLA surface 
roughness. Within the program CLA values are converted to RMS 
by multiplying by 1.25. 

Items 4 through 6 are RMS values of the slopes measured 
in degrees, of the surface asperities as measured in a traverse 
across the groove for rings, longitudinally for rollers and 
in any arbitrary direction for balls. Typical values for race- 
way and rolling element surfaces are 1 to 2 degrees. This card 
is omitted if the solution level is NPASS ■ C. 
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I FIG. 3.7 CYLINDRICAL ROLLER BEARING 

I ; FLANGE INVERSIONS 
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3. 3. 7. 2 Card Type B7B - Flange-Roller End Surface Data 

The data are identical to the data in Card Type B7A but refer 
to the flange and roller end surfaces rather than the raceway and 
roller rolling surfaces. 

Note that both card types B7A and B7B are omitted if the solu- 
tion level is NPASS = 0. 

3.3.8 Card Type B8 - Cage Data 

This card is omitted if the solution level is NPASS =0. These 
data are self-explanatory. Note that the cage weight is an input 
item. It is included mainly for future consideration of cage 
stability predictions. A nonzero value should be input to avoid 
divide checks. 

The number of degrees of freedom of the cage, MCG, is also an 
input item (either 1 or 3). When NPASS=0, cage DOF is defaulted 
to 1. See Section 2.4.1, Cage Degrees of Freedom. 

3.3.9 Card Type B9 - Shaft and Housing Fit Dimensions 

These cards are to be included only if the change in bearing 
diametral clearance with otfeiatinq conditions is to be calculated , 
i.e. if item 4 ITFIT on the Bearing Title Card is non-zero . On 
Card Type B9, tight interference fits bear a positive sign and 
loose fits, a negative sign. 

Item 3 and 6 on Card No. 9 are termed the shaft and housing 
effective widths respectively. The value specified for these 
effective widths may be as great as twice the ring width. 

Use of an effective width is an attempt to account for the 
greater radial rigidity of a shaft being longer than the ring 
that is pressed on to it, owing to the fact that the shaft deflects 
over a distance that extends beyond the ring width. In the program 
the calculated internal pressure on the ring due to its inter- 
ference fit with the shaft, is distributed over the shaft effective 
width and this (lower) pressure is used in computing the shaft 
deflection. Using double the actual width as the effective width 
is customary. 
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3.3.10 Card Type BIQ - ShaftsHousing Fit Dimensions 
These items are self explanatory. 

Note: Bearing System Components Material Properties. Card 

types Bll through B14. 

These card types define the material properties of the 
shaft, inner ring, rolling element, outer ring and 
housing, data items 1 through 5 respectively. This 
set of cards is to be included if either the bearing 
clearance change analysis is used, i.e. item 4, ITFIT, 
Bearing Title Card 2 is non zero, or if the bearing 
rings or rolling elements are not steel, i.e. item 10, 
IMT, Bearing Title Card 2 is equal to 1. If any item 
on card types Bll through B14 is left blank, the 
program inserts the appropriate value of the steel 
property. 


2.3.11 Card Type Bll - Elastic-Moduli 

This card defines the elastic modulus for the shaft, inner 
ring, rolling element, outer ring, and housing respectively. 

A default value of 204083 N/mm^ (29.6x10® PSIJ is used. 

3.3.12 Card Type B12 - Poisson’s Ratio 

This card defines the Poisson's ratio for the shaft, inner 
ring, rolling element, outer ring, and housing respectively. 

A default value of 0.30 is used. 


3.3.13 Card Type B13 - Density 

This card defines the density for the shaft, inner ring, 

rolling element, outer ^ring, and^housing respectively. 

A default value of /.806 gztmJ 15/in-^j is used. 


3.3.14 Card Type B14 - Coefficient of Thermal Expansion 

This card defines the coefficient of thermal expansion for 
the shaft, inner ring, rolling element, outer ring- and 
housing respectively, A default value of 12.24x10 ® 1/°C 
(6.8 X 10-6 i/»F) is used. 

3.3.15 Card Type B15 - Lubrication and Friction Data 

This card is omitted if the solution level is NPASS *■ 0. 
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Items 1 ind 2 


Items 1 and 2 are the amounts by which the combined thick- 
ness of the lubricant film on the rolling track and rolling 
element is increased during the time interval between the pas- 
sage of successive rolling elements, from whatever replenish- 
ment mechanisms are operative. Item 1 applies to the outer and 
Item 2 to the inner race-rolling contacts respectively. If 
Item 1 is zero or blank the mode of friction is assumed to be 
dry. 

At the present time the magnitude of the inner and outer 
raceway replenishment layers has not been correlated to lubri- 
cant flow rate, lubricant application methods and bearing size 
and speed factors. At this point then, the user is forced to 
establish proper values for the replenishment layer thickness. 

As a rough guide the following suggestions are made. 

1) To avoid starvation, the replenishment layer thick- 
nesses should be one to two times the EHD film thickness 
which develops in the rolling element raceway contacts. 

2) Because of centrifugal force effects, intuition 
suggests that the outer raceway replenishment layer 
should be several times thicker than that prescribed 
at the inner raceway. 

Item 3, XCAV, describes the percentage of the bearing 
cavitv, estimated by the user to be occupied by the lubricant. 
0<XCAV<100 . 

As with the replenishment layer thicknesses, the amount 
of free lubricant should be able to be correlated with lubri- 
cant flow rate, lubricant application methods and bearing size 
and speed factors. At this time such correlations do not exist. 
XCAV values less than five percent are recommended. 

Item 4 is the coefficient of coulomb friction applicable 
for the contact of asperities. If Item 1 and 2 are zero, then 
Item 4 serves as the coulomb friction coefficient which prevails 
in all contacts. 

Items 5 and 6 are the lubricant replenishment layer thick- 
nesses for the outer and inner ring flanges respectively. These 

appUc.ble 

for the asperity interactions at the roller end- flange contacts. 
This value should also be left blank for ball bearings. 

3.3.16 Card Type B16 

This card is omitted if NPASS title card 2 is zero or 
blank or if Item 1 cardBlBis zero or blank which implies dry 
friction. 
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This card specifies the lubricant type. If Item 1, NCODE is 
1, 2 , 3, or 4 , the Program uses preprogrammed lubricant properties 
as presented in Table 1, and no further information is required. 

NCODE Lubricant 

” 1 A specific mineral oil 

2 A MIL-L-7808G 

3 Polyphenyl-Ether 

4 A MIL-L-23699 

NCODE may also be specified as negative (-1 to -4) , in which 
case the traction characteristics of the respective lubricant 
NCODE noted above are used but the actual properties specified by 
Items 2 through 9 override the hard coded data. This option is 
most useful in specifying various mineral oils i.e. NCODE = -1. 

If items 8 and 9, AKN and ERIC are left unspecified, default 
values are set at 50.0 and 0.07, respectively. AKN and FRIC are 
only used when the NASA version of the code is exe . ted (see 
Appendix F) . 

3.4 DATA SET III - THERMAL MODEL DATA 


Appendix B has been included to aid the user in data prepara- 
tion and calculation of heat transfer coefficients required at 
input. 


3.4.1 Card Type T1 

Card type T1 is a control card. If no temperature nap is to 
be calculated, this card is to be included as a blank card fol- 
lowed by a Type T2 card for each bearing on the shaft. Card Type 
T1 contains control input for both steady state and transient 
thermal analyses. It is not intended, however, that both analyses 
be executed with the same run. 

Item l ! The highest node number (M) . The temperature nodes 
must be numbered consecutively from one (1) to the highest node 
number. The highest node number roust not exceed one hundred (100) . 

Item 2 ; Node Number of the Highest Unknown Temperature Node 
(N) . This number should equal the total number of unknown node 
temperatures. It is required that all nodes with unknown temper- 
atures be assigned the lowest node numbers. The program assumes 
that all node numbers greater than N (from N+1 to M) represent 
known boundary temperatures. 

O 

Item 3 : Common Initial Temperature (TEMP) C; The temperature 

solution iteration scheme requires a starting point, i.e., guesses 
of the equilibrium temperatures. Card Type T3 allows the user to 
input guesses of individual node temperatures. When a node is not 
given a specific initial temperature, the temperature specified 
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as Item 3 of Card Type Tl is assigned. 

Item 4 ; Punch Flag (IPUNCH); If the Punch Flag is not zero 
(0) or blank, the system equilibrium temperatures along with the 
respective node nambers will be punched according to the format 
of Card T3. This option is useful if, for instance, the user 
makes a steady state run with lubrication, and then wishes to use 
the resultant temperatures as the initiation point for a transient 
dry friction run in order to assess the consequence of lubricant 
flow termination. 
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Item 5: "Output Flaf," (lUB). If the "Output Fla?" is not 

tero’TE’e^earing program output and a temperature map will be 
printed after each call to the shaft bearing solution scheme. 

This printout will allow the user to observe the flow of the 
solution and to note the interactive effects of system tempera' 
tures and bearing heat generation rates. Two levels of bearing 
output are permitted. If lUB is 1, the rolling element output is 
not printed. If lUB is 2, full bearing output is obtained. 

Item 6: "Maximum Number of Calls to the Shaft Bearing 

Program^’’ '(ITl). ITl is the limit on the number of Thermal -Shaft * 
Bearing iterations, i.e., the external temperature equilibrium 
calculation. The user must input a non-zero integer such as 5 or 
10 in order for the Program to iterate to an equilibrium condition. 
If ITl is left blank or set to zero (0) or 1, shaft bearing 
performance will be based on the initially guessed temperatures 
of the system. The temperatures printed out will be based on the 
bearing generated heats. It is unlikely that an acceptable 
equilibrium condition will be achieved. However, the temperatures 
which result may provide better intial guesses, for a subseqmrt 
run, than those specified by the user. 

ITl also serves as a limit on the transient temperature solu- 
tion scheme, by limiting the number of times the shaft-bearing- 
solution scheme is called. Each call to the shaft-bearing scheme 
will input a new set of bearing heats to the transient tempera- 
ture scheme until a steady state condition is approached or until 
the transient solution time up limit is reached. 

Item 7: "Absolute Accuracy of Temperatures for the External 

Thermal Solution" (EPI). In the steady state thermal solution 
scheme, each calculation of system temperatures occurs after a 
call to the shaft -bearing scheme which produces bearing generated 
heats. After the system temperatures have been calculated for 
each iteration, using the internal temperature solution scheme, 
each node temperature is checked against the nodal temperature 
at the previous iteration. 

^^(K)i* ''^(N-l)i^- nodes i then equilibrium 

has been achieved and the iteration process stops. 


“ temperature of ith node at Nth thermal iteration. 
^(N-l'ii * temperature of ith node at '\*i thermal iteration. 
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I tew 8 ; "Iteration Limit for the Internal Thermal Solution" 
(IT2)*! After each call to the shaft bearing program, the 
internal temperature calculation scheme is used to determine 
the steady state equilibrium temperatures based on the calculated 
set of bearing heat generation rates. If the program is used 
to calculate the temperature distribution of a non bearing 
system it is the internal temperature scheme which is employed. 

If IT2 is left blank or set to zero, the number of internal 
iterations is limited to twenty (20), 

Item 9 ; "Accuracy for Internal Thermal Solution" (EP2) , 

The use of EP2 is explained in Section 2.1 .1 , If EP2 is left 
blank or set to zero (0), a default value of 0,001 is used. 

Item 10 ; "Starting Time" (START) is a time Tj at which 
the transient solution begins; usually set to zero (0). 

Item 11 : "Stopping Time" (STOP) is the time in s^'Conds at 
which the transient solution terminates , Tf , The transient 
solution will generate a history of the system pcrformanc-; 
which will encompass a total elapsed time of 

(T£ - Tj.) seconds 

Item 12 ; "Calculation Time Step" (STEPIN), The transient 
internal solution scheme solves the system of equations 


* tj^ ^ % LT (3.11) 

pCpV 

« STEPIN 

The user may specify STEPIN. If l®ft blank or set to zero (0), 
the Program calculates an appropriate value for STEPIN using 
the procedure described in Section 2.1.2. 

Item 13 ; "Time Interval Between Printed Temperature Maps" 
(TTIMT) seconds . The user must specify the length of time 
which will elapse between each printing of the temperature 
map. The interval will always be at least as large as the 
"calculation timestep" (STEPIN). 

Item 14 ; "Time Interval Between Calls of the Shaft Bearing 
Portion of the Program" (BTIME) . BTIME will always have a 
value larger than or equal to (STEPIN) even if the user 
inadvertently inputs a shorter inter'^al. Computational time 
savings result if BTIME is greater than STEPIN, however, 
accuracy might be lost. 
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3.4.2 Card Type T2 

Card Type T2 is reqvired, one card for each bearing if 
no thermal analysis is being performed. The temperature data 
is used within the shaft-bearing analysis portion of the pro- 
gram to fix temperature dependent properties of the lubricant 
in which cas(- the inner race, outer race, lubricant bulk 
cavity and flange temperatures are used. Note the flange 
numbering s-heme depicted in Fig. 3.5 for a cylindrical roller 
beat ii‘ For a tapered roller' bearing, the inner ring flange is considered 
flange No. i. The assembly component toiperatures at each bearing location 
are use-t in the analysis which calculates the change in bearing diametral 
clearniice from "off the shelf" to operating conditions. 


3.4.3 Card Type _T 3 

In the steady state analysis this card is used to input 
initial guesses of individual nodal temperatures for unknown 
nodes as well as the constant temperatures for known nodes, 
such as ambient air and/or an oil sump. 

In the transient analysis, Card Type T3 is used to input 
the nodal temperatures of all nodes at (START) » Tg 
at the initiation >f the transient solution. 


3.4.4 Card Type T 4 

With this card, node numbers are assigned to the components 
of each bearing, one card per bearing. With this information 
the proper system temperatures are carried into each respective 
bearing analysis. The inner race and inner ring node n;^b>rs 
may or may not be the same at the user's discretion. Similarly 
the outer race and outer ring node numbers may or may not be 
the same. 


3.4.5 Card Type t5 § T5A 

Card Type T5 is required, one card per bearing, if a thermal 
analysis is to be performed. This card designates two nodes to 
share equally each of the various types of bearing generated heat 
calculated internally by the program. For cylindrical and tapered 
roller bearings an additional card, T5A, follows immediately after 
card T5 specifying the two nodes which will equally share the heat 
generated at each of the flange contacts. For ball bearings card 
T5A must be omitted. 
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3.4.6 Card Type T 6 

This card specifies the node niunbers and the heat genera- 
tion rate for those nodes where heat is generated at a constant 
rate such as at rubbing seals or gear contacts. 

3.4.7 Card Type T 7 

This card type is used to input the nuinerical values 
of the various heat transfer coefficients which appear in the 
equations for heat transfer by conductivity, free convection, 
forced convection, radiation and fluid flow. Up to ten 
coefficients of each type may be used. Separate values of 
each type of coefficient are assigned an index number via 
card T7. In describing heat flow paths (Card Type T8 below) 
it is necessary only to list the index number by which heat 
transfers between node pairs. 

Indices 1-10 are reserved for the conduction coefficient 
X, 11-20 for the free convection parameters, 21-30 for forced 
convection, 31-40 for emissivity and 41-50 for fluid flow 
(product of specific heat, density and volume flow rate). 

As an example, for heat transfer by conduction with 
coefficient X of 53.7 watts/M°C one could prepare a card t' pe 
T7 with the digit 1 punched in column 10 and the value 53.7 
punched in the field corresponding to card columns 11-20. 

If a conduction coefficient of 46.7 were applicable for 
certain other nodes in the system one could punch an additional 
card assigning index No. 2 to the value X ■ 46.7 by punching 
a "2" in card column 10 and 46.7 anywhere within card columns 
11 - 20 . 


Rather than inputting constant forced convection coefficients, 
optionally, these coefficients can be calculated by the program 
in one of three ways. If the calculation option is exercised 
a pair of cards is used in place of a single card containing 
a fixed value of a . The contents of the pair of cards depends 
upon which of the three optional methods are used. 

Option 1) a is independent of temperature but is calculated 
as a function of the Nusselt number Nu whlth in 
turn is a function of the Reynolds number Re, 
the Prandtl number P as follows, (cf. [22] ) 
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a* NuXoil/L 

N’u » aR^pC 
e r 

where X oil is the lubricant conductivity, L is a 
characteristic length C'^’ith a unit of rocters) 
and a and b, and c are constants. 

Option 2) a is a function only of fluid dynamic viscosity 
and viscosity is temperature dependent. 

a ' cn^ 

Option 3) Q is a function of the Nusselt, Reynolds and 
Prandtl numbers -and viscosity is temperature 
dependent. 

3 ,4 . 8 Card Type T8 

This card defines the heat flow paths between pairs of 
nodes. Every node must be connected to at least one other 
node, i.e,, two or more independent node systems may not be 
solved with a single Program execution. 

The calculation of heat transfer areas is based on lengths, 

Lj and I2 input using Card Type T8. Additionally, the 
type of "surface for which the area is being calculated is indicated 
by the sign assigned to the heat transfer coefficient index. 

If the surface is cylindrical or circular the index should 
be positive, if the surface is rectangular the index should 
be input as a negative integer. 

In the case of radiation between concentric axially 
symmetric bodies, L3 is the radius of the larger body. For 
radiation between two parallel flat surfaces or for conduction 
between nodes, L3 is the distance between them. 

Fluid flow heat transfer accounts for the energy which 
the fluid transports across a node boundary. Along a fluid 
node at which convection is taking place, the temperature 
varies. The nodal temperature which is output is the average 
of the fluid temperature at the output and input boundaries. 

If the emerging temperature of the fluid is of interest, it 
is necessary to have a fluid node at the fluid outlet. At 
this auxiliary node only fluid flow heat transfer occurs and 
the fluid temperature would be constant throughout the node. 

Thus the true fluid outlet temperature will be obtained. 
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Conduction of heat through a bearing is controlled by index 
51. The actual heat transfer coefficient which contains a con- 
ductivity, area and a path length term is calculated in the bear- 
ing portion of the program. The cerm is based upon conduction 
through an average outer race and inner race rolling element con- 
tact. 


3.4.9 Card Type T9 

This card inputs data required to calculate the heat capacity 
of each node in the system. This card type is required only 
for a transient analysis and must be omirted for a steady state 
analysis . 

3.5 DATA SET IV - LOADING DATA 


If a single ball or cylindrical roller bearing is to be ana- 
lyzed, the user need not model the shaft geometry. In such a case, 
this loading card may replace all Shaft Input Data (Data Set V) . 
Data describing a dummy shaft is generated within the progreun. 

Applied loads, acting through the center of the bearing, are 
input on this card. These loads may have the form of concentrated 
radial forces (FY, FZ) , concentrated moments (MY, MZ) and a 
concentrated axial load (FX) . 

If more than one bearing is to be analyzed, this card is 
omitted, and loads are supplied with the Shaft Input Data . 

3.6 DATA SET V - SHAFT INPUT DATA 


The shaft-bearing analysis requires all lo. '*.ng to be applied 
to the shaft. The loads applied to each beari vre a product of 
the shaft-bearing solution. There is no need it the user to solve 
the statically determinate or indeterminate system for bearing 
loads. 

In the analysis the housing is assumed to be rigid. Provision 
has been allowed to input data for housing radial and angular 
spring characteristics. However, this has been done for future 
c jr.sideratlo'i of an elastic housing and is therefore currently 
unavailable. 
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The shaft input data consists of three card types: 

1) Shaft Geometry and Elastic Modulus Data 

2) Bearing Position and Mounting Error Data 

3) Shaft Load Data 
3.6.1 Card Type SI 

This card type is used to describe shaft geometry at 
up to twenty locations along tno shaft. The user must 
place his shaft in a carte .>- 5n coordinate system with the 
end of the shaft at the orig.’n and with the shaft lying 
along the X-axis. 

The shaft is described by specifying two outer and 
two inner diameters at each axial location along the 
shaft which define the diameters immediately to the left 
and the right of the X-coordinate. In this way stepwise 
or linear variations of the shaft can be handled. A 
linear variation in shaft diameter is assumed if the 
diameter to the left of one axial location is different 
than the diameter to the right of the preceding axial 
location. Complex shaft geometries may be approximated 
with a set of linear diameter variations spaced at close 
intervals . 


If an Elastic Modulus is not specified at the designated 
input location, the modulus of steel is assumed, 204083N/mm. 
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3.6.2 Card Type S2 

This card type locates the bearing inner ring on the shaft 
in the X-Y and X-Z planes. For a ball bearing, the X coordinate 
specified locates the inner ring center of curvature. For cylin- 
drical roller bearings the X coordinate locates the center of the 
inner race roller path. 

For tapered roller bearings in the strictest sense the X 
coordinate locates the point where a line from the roller center 
of gravity, intersects and is perpendicular to the inner raceway, 
with all bearing end play removed. It is sufficiently accurate 
however to allow the X coordinate to locate the center of the 
inner raceway. 

In addition to specifying bearing location, the Type 2 card 
is also used to specify housing radial and angular mounting er- 
rors. As mentioned previously, space has been reserved for in- 
putting housing radial and angular spring characteristics, how- 
ever, these characteristics are not used in the system analysis 
currently. 

Two sets of Type 2 cards may be required. The first set 
is always required and defines housing alignment errors in the 
shaft X-y plane. The second set defines the housing alignment 
errors in the shaft X-Z plane. Type 2 cards for the X-Z plane 
are required only for bearings having alignment errors in the 
X-Z plane that are different than those specified for the X-Y 
plane. The second set of Type 2 cards must be placed after the 
blank card following the Type 3 cards for the X-Y plane (see 
Appendix C) . 

The first set of Type 2 cards must contain a card for each 
bearing. The second set of Type 2 cards must give the appro- 
priate bearing number in column 10. 

3.6.3 Card Type S3 

Type 3 cards are used to specify shaft loadings at a given 
X coordinate. Loading may be applied in the X-Y and X-Z planes, 
thus requiring two distinct sets of Type 3 cards. Applied loads 
may have the form of concentrated radial forces, concentrated 
moments, linearly distributed radial forces and concentrated 
axial loads which may be eccentrically applied. If an axial 
load is eccentrically applied, the moment which results will be 
included automatically if the point of application in the radial 
plane (Y or Z coordinate) is specified in columns 71-80. Alter- 
natively, the moment generated by the axial load can be specified 
in columns 31-40. 
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Variations in distributed radial loads are handled at input 
just as shaft linear diameter variations are handled. 

Note that each set of Type 3 cards must be followed by a 
blank card. 

Also note that in order for symmetry conditions to be con- 
sidered the second Type 3 card must be void of any loading data. 

4.0 COMPUTER PROGRAM OUTPUT 


4 . 1 Introduction 


The Program Output is intended to provide the engineer or de- 
signer with a complete picture of the shaft-bearing system per- 
formance . 

In addition to the calculated output data, the input data is 
listed, thus producing a complete record of the computer run. 

Sample output of three bearing-shaft systems is included in 
Appendix D. These studies demonstrate the ability of the NASA 
version of the program to calculate performance characteristics of 
tapered and flanged cylindrical roller bearings, and to demonstrate 
the functioning of the new cage simulation algorithms and the 
single-bearinq analysis capability. The three configurations are: 

1. A system in which an input pinion is supported by tapered 
roller bearings in a straddle configuration. 

2. A system in which an input pinion gear load is supported 
by a flanged cylindrical roller bearing operating in con- 
junction with two angular contact ball bearings. 

3. A single ball bearing system operating under combined rad- 
ial and thrust loading. 

Key output items are discussed briefly below. 

4 . 2 Bearing Output 

4. 2. 1.1 Linear and Angular Deflections 

These deflections refer to the bearing inner ring relative to 
the outer ring and are defined in the inertial coordinate system 
of Figure 2.4. The bearing deflections are not necessarily equal 
to the shaft displacements since the bearing outer ring radial or 
angular mounting errors may be specified as non-zero input. 
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4. 2. 1.2 Reaction Forces and Moments 

These values reflect bearing reactions to shaft applied 
loading and outer ring mounting errors. 

When the bearing inner ring has achieved an equilibrium 
position* the summation of all bearing reaction loads should 
numerically equal the shaft applied loading. When the level 
of solution indicated by *'NPASS" ■ 2 is employed, as discussed 
in Section 5* differences between shaft applied and bearing 
reaction loads will exist but will typically be less than 101. 
This difference is a consequence of friction forces contributing 
to the reaction loads whereas the inner ring equilibrium posi- 
tion has been determined considering elastic contact forces only. 

4.2.2 Fatigue Life Data 

The Ljq fatigue life of the outer and inner raceways as 

well as the bearing are presented. The bearing life represents 
the statistical combination of the two raceway lives. These 
lives reflect the combined effects of the lubricant film thick- 
ness and material life factors. The lubricant film thickness 
life factor is described in detail in Section 3,3.1. 


4 . 2 . 2. 1 h/sigms 

The ratio h/o , also referred to as is printed for the 
most heavily loaded rolling element. The variable h, represents 
the EHD plateau film thickness with thermal and starvation 
effects considered. The variable o represents the composite 
root mean square surface roughness of the rolling element and 
the relevant raceway. 
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4,2. 2.2 Life Multipliers 

4. 2. 2. 2.1 Lubrication - This life multiplier is a function of 
h/6 at each concentrated contact- Its value’ ranges from 
0.21 for h/5 jt 0.6 to 3.0 at h/6 ^ 10. This subject is 
covered in more detail in Section 3.3.1. 


4. 2. 2. 2. 2 Material - This output simply reflects the input 
value, Again, it is covered in Section 3. 

4.2.3 Temperatures Relevant to Bearing Performance 

These temperatures fully describe the temperature con- 
ditions which affect the performance of a given bearing. If 
one of the temperature mapping options is used, the tempera- 
tures printed reflect the results of the particular option. 

If, neither temperature option was used, the list is simply a 
repeat of the input data. Note that there are separate tem- 
peratures for outer and inner raceways and flanges and ring 
temperatures. The raceway and flange temperatures are used to 
determine lubricant properties. The ring temperatures are 
used in the bearing dimension change analysis. The raceway, 
flange and ring temperatures may be the same value. 


4.2.4 Frictional Heat Generation Rate and Bearing Friction 
Torque" 


4. 2. 4.1 Frictional Heat Generation Rate 


The various sources of frictional heat generated within 
the bearings are listed. The values printed for "OUTER RACE, 
OUTER RING FLANGES, INNER R.4CE, INNER RING FLANGES, R.E.DRAG 
AND R.E. CAGE" represent the sum of the generated heats for all 
rolling elements. Additionally, the heats printed for the 
outer and inner raceways and flanges, plus the rolling clement 
cage, reflect the friction developed outside the concentrated 
contacts, i.e., the HD friction as well as the EHD friction 
developed within the concentrated contacts. The raceway and 
flange data also includes any heat generated as a consequence 
of asperity contacts when the SKP friction model is used. "R. 

E. DRAG" should be interpreted as the heat resultinq from lubri- 
cant churning as the rolling elements plow through the air-oil 
mixture. 
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4. 2, 4. 2 Torque 

The torque value is calculated as a function of the total 
generated heat and the sum of the inner and outer ring rotational 
speeds. The intent is to present a realistic value of the torque 
required to drive the bearing. Under .conditions of inner ring 
rotation the torque value reflects the torque required to drive 
the inner ring . 


4 . 2 . S EHD Film and Heat Transter Data 

4. 2. 5. 1 EHD Film Thickness 

These values refer to the calculated EHD plateau film 
thickness at both contacts of the most heavily loaded rolling 
element and include the effects of the thermal and starvation 
reduction factors. 

4 . 2 . 5 . 2 Starvation Reduction Factor 

These factors give for the inner and outer ring contacts, 
the reduction in EHD film thickness ascribable to lubricant 
film starvation according to thv'i methods of Chiu, [11]. 

These factors pertain to the EHD film thickness for both 
the inner and outer race contacts of the most heavily loaded 
rolling elements, but are applied to the respective inner and 
outer race film thickness for each rolling element in the bear- 
ing. 

4.2. 5.3 Thermal Reduction Factor 

These factors are calculated according to the methods of 
Cheng, [lO] and pertain to the EHD film thickness for both 
the inner and outer race contacts of the most heavily loaded 
rolling elements, but are applied to the respective inner and 
outer race film thickness for each rolling element in the bear- 
ing. 


4. 2. 5.4 Meniscus Distance 

These values are calculated according to the methods of 
Chiu, [11] and pertain to the EHD film thickness for both 
tho inner and outer race contacts of the most heavily loaded 
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rolling elements, but are applied to the respective inner and 
outer race film thickness for each rolling element in the bear- 
ing. 

4. 2. 5. 5 Raceway- Rolling Element Conductivity 

These data reflect the amount of heat transfer between 
rolling element and raceway for each degree centigrade differ- 
ence between the two components. These data reflect the average 
of all outer and inner contacts respectively. (See Sections 2.1.4 

2. 1.4.1 for a discussion of the calculation procedure.) 

4.2.6 Fit and Dimensional Change Data 

4. 2.6.1 Fit Pressures 

These data refer to the pressures built up as a conse- 
quence of interference fits between shaft and inner ring and 
housing and outer ring. Pressures are presented both for the 
standard cold-static condition (160C) and at operating 
conditions . 

4. 2.6. 2 Speed Giving Zero Fit Pressure (B etween the shaft and 

inner ring ) ~ 

This is a calculated value based upon operating conditions 
and provides a measure of the adequacy of the initial shaft 
fit. 

4. 2. 6. 3 Clearances 

"Original" refers to cold unmounted clearance which is 
specified at input if the diametral clearance change analysis 
is executed. "Change" refers to the change in diametral 
clearance at operating conditions relative to the cold unmounted 
condition. A minus sign indicates a decrease in clearance. 
"Operating" refers to the clearance at operating conditions. 

For all types of ball bearings the decrease in clearance can 
be combined with the initial diametral clearance, and the free 
operating contact angle at operating conditions may be 
calculated. Note that the change in clearance should be com- 
pared against the diametral play of the split inner ring ball 
bearing in order to determine if the possibility for three 
point contact exists. The Program does not account for three 
point contacts even though the change in clearance mignt suggest 
that three point contact is obtained. 
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4.2.7 Lubricant Temperatures and Physical Properties 

The lubricant properties, particularly the dynamic vis* 
cosity and to a lesser degree the pressure viscosity coefficient, 
are heavily temperature dependent. These factors enter the EHD 
film thickness calculation and the HD and EHD friction models. 

The lubricant is assumed to be at the same temperature as the 
relevant raceway. As noted elsewhere, these temperatures may 
be either input directly or calculated by the Program. 

The physical properties printed are seif explanatory. 

The units are enumerated. 

4.2.8 Cage Data 

4 . 2 . 8 . 1 Cage-Land Interface 

The cage data indicate the performance parameters at the 
interface between the cage rail and the ring land on which the 
cage is guided. The torque, heat rate and separating force 
require no explanation. The eccentricity ratio defines the 
degree to which the cage approaches the ring on which it is 
guided at the point of nearest approach. The radial displace- 
ment cf the cage relative to the bearing axis is divided by one 
half tne cage- land diametral clearance. An eccentricity ratio 
of one indicates cage-land, contact. A ratio of zero indicates 
that the cage rotation is concentric with the bearing axis. 

Only the cage-land and rolling element pocket forces are 
considered in determining the cage eccentricity. The cage 
weight and centrifugal force which result from the eccentricity 
although available, are not considered in the analysis. The 
omission of these considerations helps reduce convergence 
problems, 

4. 2. 8. 2 Cage Speed Data 

Cage speed data present the comparison between t a cage 
speed calculated based upon the quasidynamic equilibrium con- 
siderations and the speeds calculated with raceway control 
theory for ball bearings and tj\e epicyclic speeds of the 
roller bearing components. 

4.3 Rolling E lement Output 

4.3.1 Rolling Element Kinematics 
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4 . 3 . 1 . 1 Rolling Element Speeds 

All of the rolling element speeds tend to vary from posi- 
tion to position when the bearing is subjected to combined 
loading. 

The total rolling element speed is with reference to the 
cage and represents the vector sum of the three orthogonal 
components . 

4. 3.1.2 Speed Vector Angles 

The rolling element speed vector angles, Arctan 
and Arctan (wj/wx) are presented in order to show a clearer 
picture of the predicted ball kinematics. The ball speed 
vector tends to become parallel with the bearing X axis with 
increasing shaft speed and decreasing contact friction. 

4.3.2 Rolling Element Raceway Loading 

4. 3.2.1 Normal Forces 

The normal forces acting on each rolling element are 
printed. The rolling element race normal forces are self 
explanatory. The cage force is calculated only when the 
friction solution is employed and is always directed along the 
rolling element 2 axis. If the rolling element orbital speed 
is positive, a positive cage force indicates that the cage is 
pushing the rolling element, tending to accelerate it. Cage 
force is a function of rolling element position within the 
cage pocket. Its magnitude is derived using hydrodynamic lubri- 
cation assumptions, when the distance between the rolling element 
and cage web is large, and EHD assumptions when the separation 
is of the order of the EHD film thickness or when rolling 
element- cage web interference exists fSKF version). The NASA 
version (see Appendix F) solves for cage forces using hydrodynamic 
lubrication assumptions only. 

4 . 3 . 2 . 2 Hertz Stress 

The stress printed represents the maximum normal stress at 
the center of each ball race contact or at the most heavily loaded 
slice of the roller raceway contact. 
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4. 3. 2. 3 Load Ratio 

asp -tot- 

If the EHD film thickness is small compared to the RMS compo- 
site rolling element-race surface roughness, the rolling element- 
race normal load will be shared by the EHD film and asperity con- 
tacts. The load ratio reflects the portion of the total load 
carried by the asperities. 

When the NASA version of the code is executed (see Appendix F) , 

^asp^®tot=0, since asperities are not accounted for in the Allen 
Traction model (71 . 
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4. 3. 2. 4 Contact Angles 

A ball bearing, subject to axial loading, misalignment 
or mounted such that the inner ring is always displaced axially 
relative to the outer rings, (i.e. a duplex set of angular 
contact ball bearings) will have non-zero contact angles. At 
low ball orbital speeds the inner and outer race angles will be 
substantially the same. At high speeds, ball centrifugal force 
will cause the o iter race contact angle to be less than the 
inner race angle. 


4.3.3 Roller End-Flange Contact Data 

For a taperei roller bearing a single set of roller 
end-flange data is printed. For the cylindrical roller bearing, 
which may have jp to four flanges, the program examines the 
data and prints the results for the most heavily loaded outer 
and inn r ring flanges. In addition to the data listed belcw 
for a cylindrical roller, the semimajor contact axis as well as 
the concentrated contact and hydrodynamic heat generation rates 
are printed for the tapered roller. 

4. 3.3.1 Normal Fo rce 

The interference between the roller end and flange is 
dcf^rminei from che solution for the relative ?ring and rolling 
element positions. Hertz theory is used to calculate the load 
which results from this interference. 

4. 3. 3. 2 Hertz Stress 


The Hertz stress printed is the maximum normal stress 
which occurs in the contact. 

4. 3. 3. 3 EKD Film Thickness 


The plateau film thickness is calculated using the 
Archard-Cowking [3] equation for point contacts or the Dowson- 
Kigginson[4] equation for line contacts, when the "SKF" version of 
the )rograin is executed (see Appendix F) . The ’’NACA” version uses the 
Loewenthal Pjdel [6J to calculate film thickness. Either result is 
then modified to account for starvation [12] and thermal [11] reduction 
effects. This modified film thickness is printed. 

4 . 3 . 3 . 4 Sliding Velocity 

Til'- sliding velocity is defined as the difference between the 
flange and roller end linear velocities at < he center of the contact. 
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4. 3, 3. 5 Rolling Velocity 

Rolling velocity is defined as one half the sun of the 
roller end and flange linear velocities at the center of the 
contact. 


4. 3. 3. 6 Contact Ellipse Semiminor Axis 

To help assess the severity of the roller end flange 
contact and the possibility for edge loading* the semiminor 
axis of the contact is printed, 

4.4 Thermal Data 


As in the case for bearing output* all o.f the input data is 
printed. The calculated output data is presented in the form 
of a temperature map in which a node number and the respective 
node temperature appear. The appearance of the steady state 
and transient temperature maps are identical. The transient 
temperature map also includes the time (T) at which the 
temperature calculations were made. 

4. 5 Shaft Data 

These data simply reflect the input information. 

4.6 Program Error Messages 

4.6.1 From Subroutine ALLT 

"Steady State Solution with (EPl) degrees accuracy was not 
obtain - after (ITl) Iterations". 

This message pertains to the external temperature itera- 
tion scheme in which system temperatures and bearing generated 
heats are being solved for an equilibrium condition. 

4.6.2 From Subroutine SHADE 

"It was not possible to obtain the change of clearance 
with an accuracy of (ERFITJ times the rolling element diameter 
(JTFIT) iterations": 
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This message pertains to the bearing diametral clearance 
change iteration scheme. The solution may be converging in 
which case the number of iterations (ITFIT) should be increased. 
This can be checked with an NPRINT ■ 1 intermediate printout. 

The intermediate rint may indicate that the solution is 
oscillating. The most likely cause of oscillation is the 
alternate predict! m of bearing preload with all rolling elements 
loaded, and then in the next iteration, only a subset of the 
rolling elements loaded. This problem can usually be over- 
come by either of two methods. 

1) In subroutine FIT remove the GO TO 20 statement. 

This will cause the inner ring load distribution 

to have no effect on the change in diametral clear- 
ance. 

2) The solution can be damped by redefining the solution 
damping factor FA, such that it would take on a value 
0<FA^. FA is presently set to 1 in subroutine SHABE. 

If this damping technique is used, the number of FIT 
iterations should be increased as the value of FA is 
decreased. An upper limit of 10 iterations is recom- 
mended. 

4.6.3 From Subroutine SOLVXX 

1) "SINGULAR SET OF EQUATIONS" 

This message night occur when the thermal input data 
is not input properly. 

2) "THE LIMIT FOR NUMBER OF ITERATIONS IS REACHED" 

This message might occur either during a steady state 
temperature solution or bearing solution. Before 
increasing the number of iterations check the equa- 
tion residue values. If they are low, the solution 
may be good enough. 

3) "THIS IS THE BEST WE CAN DO. IT MAY BE USEABLE" 

This message reflects the fact that the next iteration 
will result in divergence. The iteration procedure is 
thus terminated. The equation residue values should 
be checked, if they are low, the solution may be useable. 
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As noted above the occurrence of messages 2) and 3), do 
noC necess irily mean that the solution is not good. Generally 
the messages indicate that the solution has not converged as 
tightly as the user has requested. 

Note: The XX suffix on SOLV specifies the version of SOLV 
contained in the user's program. As of this writing the current 
version is SOLVl 3. The suffix is changed each time improve- 
ments are made which require a change in the SOLVXX calling 
sequence. 

4.6.4 From Subroutine INTFIT 

“Sjpgular matrix on tight shaft fit" 

"Singular matrix on loose shaft fit" 

These messages reflect an erro^ in the input data usually 
as a consequence of inconsistent component diameters » such as 
the shaft inside diameter being greater than the outside 
diameter. 
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5.0 GUIDES TO PROGRAM USE 

The Coaiputer Program is a tool. As with any tool the 
results obtained are at least partially dependent upon the 
skill of the user. Certainly the economics of the Program 
usage are highly dependent upon the user's technical need 
and discriminate use of Program options. 

Some general guides for efficient use of the Program are 
listed below: 

1. Attempt to use the lowest level of solution possible . 
For instance if the prime object of a given run is to 
obtain bearing fatigue lives, execute only the elastic 
solution (NPASS * 0). If an estimate of bearing 
frictional heat is required, execute the low level 
friction evaluation (NPASS » 1). Execute the friction 
solution (NPASS * 2) only if rolling element and cage 
kinematics are of interest. Execute the highest level 
of solution (NPASS - 3) if kinematics are of interest 
and the bearing reaction loads deviate substantially 
from the shaft applied loading, i.e. a deviation 
greater than ten percent, 

2. Attempt to input bearing operating diametral clearance 
rather than calculate it. Or, execute the diametral 
clearance change analysis once for a group of similar 
runs and use the output from the first run as input 

to the subsequent runs omitting the clearance change 
analysis . 

3. Attempt to input accurate operating temperatures 
rather than calculate them. 

4. The more non-linear the problem the more com puter time 
required to solve it. In the bearing friction solution 
large coefficients of friction seem to increase the 
degree of nonlinearity. In the thermal solttions, if 
possible, eliminate nonllnearities by omitting 
radiation terms and by using constant rather than 
temperature dependent free and forced convection 
coefficients. 

5. In the transient thermal solution, space the calls to 
the shaft-bearing solution TBTIMeI to as large an 
i nterval as prudently possible . Be careful however , 

Too long and interval will produce large errors in heat 
rate predictions. 
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6. In the steady state thermal analysis attempt to esti - 
mate nodal temperatures on a node by node pa^is . 

Nodes which are heat sources shoula have higher 
temperatures than the surrounding nodes. 

The above suggestions are intended to encourage the use 
of the Program on a cost effective basis. The intent is not 
to discourage the use of important program capabilities but 
to emphasize how the program should be most effectively used. 

It is suggested that the user take a simple, axially 
loaded ball bearing problem and execute the program through 
the full range of options beginning with a frictionless solu* 
tion proceeding to the three levels of friction solution with 
a low (0.01) and high (0.1) friction coefficient. The dia- 
metral clearance change analysis and the thermal solutions 
should also be executed on an experimental basis. This exercise 
will provide the user with some insight into economics of the 
Program usage on his computer as well as the results obtained 
from various levels oi solution of the same problem. 

It is also suggested that a constant user of the program 
should study the hierarchical Program flow chart Appendix A, 
along with the Program listing to gain an appreciation of. the 
program complexity and the flow of the problem solution. The 
Program is comprised of many small functional subroutines. 
Knowledge of these small elements may allow the user to more 
ea.sily piece together the philosophy of the total problem 
solution. 
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Flow Chart 


The hierarchical flow chart presents the program structure, 
listing the program elements in the order in which they would 
be called to solve the shaft-bearing dynamic, as well as steady 
state and transient temperature distribution problems. The 
various solution loops are indicated, as well as notes which 
indicate the functions of various subroutine groupings. 

Each line in the flow chart represents a program element, 
subroutine, function or the main program ALWAYS. The call of 
one subroutine by another is denoted by indenting the called 
subroutine relative to the routine doing the calling. As an 
example, subroutine SKF calls subroutines FLAGS, TYPE, PROPST, 
LUPROP, LUBCON, DATOT, CNVRT, CONS and SPRING. Subroutine 
CONS calls CONBRl, BCON, TCON, CRCON and C0NBR2. BCON calls 
ABDEL. Both TCON and CRCON call ABDEL and SLICES. 

The first mention of a subroutine within the flow chart 
includes the entire list of subordinate program elements. At 
subsequent calls to that subroutine the list of subordinate 
elements is omitted. As an example the first call to subroutine 
AXLBOJ is followed by the subordinate elements, JMVIKT, SNITMT, 
NUMLOS, DUBSIM, MEIE, MEIL and SIMQ. After the call of AXLBOJ 
from INDEL, the subordinate elements are not listed but are, 
nevertheless, employed. The list of subordinate program ele- 
ments are omitted in repeated calls of subordinate GUESS, BEAR, 
SOLVXX and DELIV3 as well as AXLBOJ. 

As noted earlier, rolling equilibrium is calculated, first 
without, then if required, with friction forces included. 

Whether or not friction is considered is highlighted with the 
words Frictionless or Friction. 

If the Program is too large to fit in its entirety on the 
user’s computer, segments of the program may be ’’overlaid”. 

For this purpose the Program is subdivided into ten (10) 
modules which can be sequentially ’’overlaid”. 

The Program segments SKF, TEMPIN, and GUESS all per- 
form initiation functions and with the exception of GUESS, are 
called only once per program execution. 

The real problem solving portion of the program is e’- jdied 
in se^ent ALLT. Within this segment the shaft bearir ^ solu- 
tion is obtained th-^^ugh the call to SHABE, then the steady 
state or transient .eraperature distributions are obtained. 
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This scheme is repeated until the end objective, steady state 
thermal equilibrium or time up for the transient scheme, is 
realized. 

The nonlinear equation solver S0LV13 is central to the 
program and deserves special discussion as related to the flow 
chart. The first call to S0LV13 is BEAR. Only for this 

first call are all of the S0LV*13 subc.ainate subroutines 
listed as noted earlier. These include INSOLV, EQS, PARDER, 

SIMQ, EQCHEK, ERWRIT and ERCHEK. In the subsequent call to 
S0LV13 in which the steady state temperatures are being calcu- 
lated, the above listed subroutines are again called but 
these calls with the exception of EQS are not listed on the 
flow chart. 

E(^S is the name given by S0LV13 to a subroutine which sets 
up the system of equations to be solved. EQS is brought into 
S0LV13 through the argument list. When the bearing equations 
are being solved, subroutine BRGGEQ is brought into S0LV13 
and within S0LV13 is referenced by the name EQS. When the heat 
transfer equations are being solved as a consequence of the 
call of S0LV13 from ALLT, NET is brought into SOLV13 and is 
referenced as EQS. 

Since storage and execution costs rise dramatically when the 
program is run at levels higher than NPASS=1, it is recommended that 
the program be executed at the lowest solution level possible. If 
a level 2 or 3 shaft-bearing solution is desired for bearings with 
more than 16 rolling elements, a new executable version must be cre- 
ated by expanding the blank common used in' subroutine BEAR. Increas- 
ing the length of the C array to 33489 elements enab' the program 
to simulate bearings with up to 30 rolling elementf 
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VERSION) 


JViad and lvI bonring crK^ bearing solution 
oontr I dsU 


Stt mtftrlsl and lubricant propartias and 
calculalu oorciUaiLo !br' UM|>craLuro dcix^ndcncy 
calculations 

Writ# bearing input and bardcodad, preset datr 


Caloulate bearirwi related oonatanta 


Paad and write theroal and thernal solution 
ocntrol dau and oaloulaU heat transfer 
ooefficients 


l^ed and write shaft geoectry, lending and 
bearing poaition data 


Calculate shaft defleotior eonstanta 


Calculate shaft influence coefficients 


Quess vtlues of noiung element and cag^' 
variables, posi . ons and speeds 


I 


Bogin tlie solution of the steady state transient 
teoperature dependent shaft-bearing 

Calculate bearir<4 dianetral clearance 

Establish iteration schsas to satisfy shaft 
and tearing inner ring equillbriua 
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Of WORQUAUrt 
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PREPAY 
I INITX 
I i UMLOAD 

I k UNIX)DT 

I KAXMIN 
SOLVli 
I INSOtv 
I rxjaiiic 
I I DAPiPCO 

I I ERCHEK 

I E3QS « BnCGLXJ (WITHOUT FROTION) 
i I BQBAa (BAa BEARING) 

1 I f BCCTRL 

t I I BALLXN 

I ) I BALLEQ 

I { EX7TAPR (TAPERn) ROUINC BEARING) 

I i I ucxrniOL 

i I I TAPIN 

I I I TAPED 

I I t I TPNORM 

j ! I I FTHORH 

I I I I PUSH 

t I EQCYL (Cini^ *:aL ROIXER BEARL 4G 
I I I BGCr<L 

I 1 I ROUJN 

1(1* FLNDEF 
ill XQ 
I I I 1 TPMOfW 

I I I I rMom 

I I I I I 

1 I I I I PULL 

I PARDER 

I I EOS ■ BRGGEQ 

I SIMQ 

I BQCHE3C 

I I EQS > BRGGEQ 


EiUbllah Itaratlon schane to misfy rollir 4 
•Ictient equilibrium. Solve rolling elumunt 
•quili*^ Ijn, bearir\s Ly t»arlng, one eleiiisnt 

at a Ut a 


Evaluat'^i ball-f*aceway contact loads, ball 
centriiugal force and ball equilibria equations 


Evnluito tnporod rollorTnccwy ond flnnpT^ 
contact loads, centrifugal force and gyroooment 
and roller equlllbrlue equations 


) Evaluate cylindrical roller-raceway and flange 
contact loads, centrifugal force and roller 
equllibrluB equations 


Calculate partial Jerlvatlves of rvllirg elenant 
equilibrium equations with respect to rolling 
element positions and calculate new positions 


r BFILL 
SHAPA 
) FIUB 
SIMEQ 


P Vise ? 

I ALPHAO 
DRACNO 
STCCN 


EVALUT (If NPASS s 1) SKIP II' NPASS « 2 OR 3 
I PREPAR 

I BHXEQ (WITH FRICTION) 


EQBALL 
1 BGCTRL 
I BALLtN 
I BALLEQ 


FMIX/NASA 
I TINT 
I FILM 
I THERTl 
I STARFC 
I HDFRIl 
I 1 FBAR 

I FRINT/NASA 
I ALLEN 

I ' IRICTN 


CACESP 
CAGEEQ 
I CGLAND 
I i CCDRY 

I OCMT 
I CCiiE 

I I CGBALL 

I I CGROLL 

I SUMCGL 
I SUMFE 


I i 1 I r 

I I I I 1 


Sm rolling element-li..ier ring forces and ocments 


Calculate bear)* r fatigue life 
Add partial derivatives of bearing inner ring 
force with respect to inner ring displacement 
to the shaft equilibrium equation end obtain 
new inne! ring positions 

After shaft equilibrium is satisfied, caieuintc 
bearing reaction forces at the bearlr^ equilibrium 
posiUona and astlmate a r>^ set of rolling element 
spems 

E^iuate temperature dependent lubricant 
prope'"tles and constants. 


Evaluate bearing performance with estimated 
rolling element speeds. 


Evaluate ball -raceway contact loads. 


Evaluate ball -raceway EH) film tMcur.ess, 
inlet aril concentrated contact friction 
forces and heat generation rate. 


Evaluate caga speed, Ul-cage normal mvi 
friction forces and heat ratea plus cage 
rail -ring land normal and friction forces 
plus heat genaration rates. 
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I nnnAX 
I 1 OHAX 


! HurnAc 

I 1 HUT 

EX3TAPR 
I BGCTRL 
I TAPIN 
1 TAPBQ 
I 1 TPNORM 

I I FLNORM 

I 1 I PUSH 

1 1 I pua 


fmhr/nasa 

I SPEEDS 
i riLM 
I THERFl 
I STARPC 
I HRFRIC 
I SPEEDS 
i ALLEN 
1 FRICTN 

FLMIX /NAwA 
1 FILM 
) THERFl 
\ STARFC 
I HDFRIl 
t I FEAR 
I ALLNPT 
I FRICTN 


I EQCYL 
I I BGCTRL 

I I ROLLIN 

111 FLHDEF 
I I RaLEQ 

III IPNOPM 
II; FLNORM 
I I I i PUSH 
111 FMIXR /NASA 
III RMDC/NASA 

1 I 

I I CAGSPD 

I 1 CAGEEQ 

! I I BRGAX 

I I I LIFE 

I I HUTRAC 

I I I HUTFLN (IF NPASS = 1, RETURN TO AUT.) 

BEAR (IF NPASS = 2) NPAbi = 3 

I PREPAR 

4 1 savi3 

l I EOS 5 BRGGEQ (WITH FRICTION) 

SONRl (IF NPASS = 3) 

I BEARC 
3 I I BEAR 

I 1 I PREPAR 

I I#** I SOLV13 

111 1 EOS 5 BRGGEQ (WITH FRICTION) 

I I I SUMF 

I I I I SIHC 
1111 SUMFL 


DELIV3 /SKF 
1 TITLE 
I RITE 
I RITEl 
I REOUT3 
I i RITE2 

I REFLOT 
1 I RITE4 

I RITE2 

I FLNGOT 


Evaluate ball, position depondont, aocoloratlona, 


Calculate bearing fatigue life 


C&loilate ballHracoway heit transfer coefficients 


Evaluate roIler-racMay and roller-flanm contact 
loads 


Evaluate roller raceway EHO film thickness inlet 
and ooncentrated contact fk*lction forces and 
heat generation rate 


Evaluate roller-flange EKD film thickness inlet 
and concentrated contact friction forces and 
heat generation rate 


Evaluate tapered roller bearing cage speed, cage 
forces, roller accelerations, bearing fatigue 
life, roller-raceway and roller-flange heat 
transfer coefficients 


Evaluate cylindrical roller-raceway and flange 
normal and friction forces and heat generation 
rates plus EHD film thicknesses 


Evaluate cylindrical roller bearing cage speed, 
cage forces, roller acceleration, beai'ing fatigue 
life, plus roller-raceway and roller-flange heat 
transfer coefficients 


Solve tr^ rolling element and cage equilibriijn 
equations with friction included, using the 
inner ring positions established by t*he 
frictionless solution 

Solve shaft, rolling element and cage equilibrlun 
equations with friction included 


Write shaft-bearing output 


origmwu. page is 
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FILLGT 

sa.vi3 

I BQS = NET (FOH STEADY STATE TUEJWAL SOLUTION) 
1 I NETEET 
DELIV3 
■IMAF 


STEPMA 

I NET (FOB THE TRANSIENT TOEfMAL SOLUTION) 
1 1 NETEET 

DELIV3 
NET 

I NETFET 
TNAP 

[ Daiv3 

TMAP 


Solve for the sLcady state temperature 
distribution. Write shaft -bearing output 
arKl temperature mr»r> if dr:;iix>d. 


Determine an appropriate time step and 
solve for the transient tenperature 
distribution. Write shaft -bearing output 
and tempereture map if desirsd. 


Write output for the final solution 


NPASS s 0 SHAFT At4D BEARING INNER RING EQUILIBRIljM ARE 

S/iTISFIrD CONSIDERING ELASTIC CONTACT FORCES. 

NO LUBRiaiTOK OR FRICTION EFFECTS ARE 
CONSIDERED. 

NPASS « 1 SHAFT* AND BEARING INNER RING EQUILIBRIUM ARE 

SATISFIED CONSIDERING ELASTIC CaTTACT FORCES. 
i.iniRTrATTON AND niTrTUiN iJ-r-Tciv. AiJi*. crHY:uni:nF3) 
USING RACEWAY CONTROL (BALL BEARINGS) O.R 
EPICYaiC (ROLLER BEARINGS) ASSUMPTIONS TO 
ESTIMATE ROLLING ELEMENT AND CAGE SPEEDS. 

NPASS a 2 SHAFT AND BEARING INNER RING EQUILIBRIUM ARE 

SATISFIED CONSIDERING ELASTIC CONTACT FORCES. 
USING THE INNER RING POSITIONS THUS OBTAINED. 
ROLUNG ELEMENT AND CAGE EQUILIBRIUM ARE 
DETERMINED CONSIDERING FRICTION. 

NPASS « 3 COFSPLETE SOC.UTTON. SHAFT AND BEARING INNEH^ 

RING PLUS HOLING ELEMENT AND CAGE EQUILIBRIUM 
ARE DETERMINED CONSIDERING ALL ELASTIC AND 
FRiaiON PCfflCES. 


A;7 



AT81D040 


APPENDIX B 


HEAT TRANSFER COMPUTATION NOTES 
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B.l BACKGROUND 

The temperature portion of Program AT81Y003 is designed 
to produce temperature maps for an axisymmetric mechanical system 
of any geometrical shape. The mechanical system is first approxi- 
mated by an equivalent system comprising a number of elements of 
simple geometries. Each element is then represented by a node 
point having either a known or an unknown temperature. The environ- 
ment surrounding the system is also represented by one or more 
nodes. With the node points properly selected, the heat balance 
equations can be set up accordingly for the nodes of unknown 
temperature. These equations become non-linear when there is 
convection and/or radiation between two or more of the node points 
considered. The problem is therefore reduced to solving a set 
of linear and/or non-linear equations for the same number of 
unknown nodal temperatures. It is obvious that the success of 
the approach depends largely on the physical subdivision of 
the system. If the subdivision is too fine, there will be a large 
number of equations to be solved; on the other hand, if the 
subdivision is too crude, the results may not be reliable. 

In a system consisting of rolling bearings, for the sake 
of simplicity, the elements considered are usually axially 
symmetrical, e.g., each of the bearing rings can be taken as an 
element of uniform temperature. For an element which is not 
axially symmetrical, its temperature is also assumed to be 
uniform and its presence is assumed not to distort the uniformity 
in temperature of a neighboring element which is axially symmetri- 
cal. That is, the non- symmetrical element is represented by 
an equivalent axially symmetrical element with approximately 
the same surface area and material volume. This kind of approxi- 
mation may seem to be somewhat unrealistic, but with properly 
devised equivalent systems, it can be used to solve complicated 
problems with results satisfying some of the important engineer- 
ing requirements. 

The computer program can solve the heat-balance equations 
for either the steady state or the transient state conditions 
and produce temperature maps for the mechanical system when the 
input data are properly prepared. 

b.2 basic equations 

B. 2 . 1 Heat Conduction 

The rate of heat flow q^. 4 (^Vf) that is conducted from node 
i to node j may be expressed qy. 
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and tj are the temperatures at i and j, respectively, j 
the area normal to the heat flow, lij the distance * 

(m) and the thermal conductivity oetween i and j, (W/m°C). 

Assuming that the structure between point i and j is 
composed of different materials, an equivalent heat conductivity 
may be cal ulated as follows: 




J 


Fig. B-1 

A.. . Vl * ^2 *2 


Aj. • Aj . Aj 


Fig. B-2 


X . . « 
ij 




VS * S/S 


f . . ■ A, + 1. 

ij 1 2 


The calculation of the areas will be discussed in Section B.2.5. 
B.2.2 convection 

The rate of heat flow that is transferred between a solid 
structure and air by free convection may be expressed by 
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where 


in which 




a- .• A. . It. - t. 


1.25 


vi.j “i.j ''i.j '^i 


• SIGN (t. - t.) 


SIGN 


. f 1, if - tj) > 

1-1. ifCti - tp< 


2.5 • 10’^ W/m^- (degC)^*^^ for hot surfaces facing upward 
and cold surfaces facing downward 

‘ij ■V .4 • lO"*' W/m^- (degC)^'*^ for hot surfaces facing downward 

and cold surfaces facing upward 

[1.8 • lO'^ W/m^- (degCj^**^ for vertical surfaces 

For other special conditions, o^j must be estimated by referring 
to heat transfer literature. 

The rate of heat flow that is transferred between a solid 
structure and a fluid by forced convection may be expressed by 


■ “i.j *i.j 


in which is the heat transfer coefficient. 

Now, with ®ij » introduce the Nusselt number 


the Reynolds number 
and the Prandtl number 


where 


aL 

T 


ML 

V 

p V 

£ 

X 


L is a characteristic length which is equal to the diameter 
in the case of a cylindrical surface and is equal to the 
plate length in case of a flat surface (m) 
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U is a characteristic velocity which is equal to the difference 
between the fluid velocity at some distance from the surface 
and the surface velocity (in/sec ) 

A is the fluid thermal conductivity (W/M° O 

V is the fluid kinematic viscosity (M^/sec) 

p is the fluid density (kg/m^) 

Cp is the fluid specific heat(J/kg°C) 

For given values of R_ and Pj- the Nusselt number and 
thus the heat transfer coefficient may be estimated from one 
of the following expressions: 


Laminar flow along a flat plate: Rg < 2300 

- 0.323 /R^ ♦ 

Laminar flow of a liquid in a pipe: 

^ 

where D is the pipe diameter and L the pipe length 


Turbulent flow of a liquid in a pipe: 
» 0.027 . R°-® . \ 

Gas flow inside and outside a tube: 


\ = 0.3 R°*^^ 


Liquid flow outside a tube: 


N 


u 


0.6 R 


0.5 


p0.31 

r 


Forced free convection from the outer surface of a 
rotating shaft 

Ny « 0.11 [o.5 Rg^ . 

where the Reynolds number Rg is developed by the shaft 
rotation. 

„ OUT 

Re “ -V 
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in which w is the angular velocity (rev/sec) 

D is the shaft diameter (m) 

The average coefficient of forced convection to the lubri- 
cating oil within a rolling contact bearing nay be approximated 
by. 


a - 0.0986|’^[l ^ 

using + for outer ring rotation 
- for inner ring rotation 

in which N is the bearing operating speed (rpm) 

D is the diameter of the rolling elements (mm) 
dm is the bearing pitch diameter (mm) 

B is the bearing contact angle (degrees) 


B.2.3 FLUID FLOW 


The rate of heat flow that is transferred from fluid 
node i to fluid node j by fluid flow is 








Vjj is the volume rate of flow from i to j . It must be observed 
that the continuity of mass requires the following equation to 
be satisfied 


provided the fluid density is constant. The summation should 
be extended over all nodes i within the fluid which have heat 
exchange with node j by fluid flow. 

B.2.4 HEAT RADIATION 

The rate of heat flow that is radiated to node j from 
node i is expressed by 

■’RiJ ■ ‘l.j (CV273)^-Ctj*273)*) 

where 

T^ • temperature of node i in °C. 

T^ ■ temperature of node j in ®C. 
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and the value of the coefficient <5, j depends on the geometry 
and the emissivity or the absorptivity of the bodies involved. 

For radiation between large, parallel and adjacent surfaces 
of equal area, Aj j and emissivity, j, j is obtained from 
the equation * 


where 


0 

I 


1.1 1,1 



the Stefan-Boltimann constant, 


is 


5.76 


10*® W/m^ /(degK)^ 


For radiation between concentric spheres and coaxial cylinders 
of equal emissivity, is given by the equation 


1, 1 


. oA 


ijJ_ 


11 


* * tl- *i i> A^. 

i.i 


where o is as above, A^ j is the area of the enclosed body 
and A*i^j is the area oi the surrounding body, i.e. Ai^j 

Expressions for "Sijj that are valid for more complicated 
geometries or for different emissivities may be found in the 
heat transfer literature. 


B.2.5 CALCULATION OF AREAS 


In the case of heat conduction heat transfer in the axial directio 
IS given by the equation (Fig. B-3). 


^i.i 


2irr.„ • 
m 


Referring to the input instructions, Section 3.4.7, but recalling 
L must be input in mm not m. 


4 • ’^m 


^1 " ^2 


L^ - Ar 


^2 ’ "l 
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In the case of heat transfer in the radial direction, j 
is obtained from the expression * 



2TT r^ . H; L, - r^; L, - H 
m X m £ 


and similarly for the radiation term above 


A*. . - 2ir r* H 
1,1 m 


L 


3 


r 


* 

m 


L 


2 


H 


in which H is the length of the cylindrical surface; where heat 
is conducted between i and j , ri^i is given by the same equation 
as above (Fig. B-4 (a) ); where heat is convected between'i and j. 



Fig. B-3 


Fig. B-4(a) 
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B . 3 . 1 TRANS I ENT ANALYS I S 

For the transient analysis all of the data pertaining to 
the node to node heat transfer coefficients must be provided by 
the input. Additionally, the volume and the specific heat at 
each node is required. 
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SKK COMPUTER PROGRAM SHABERTH 
SAMPLE OUTPUT 


The SHABERTH output samples which are displayed on Pages D:3 
to D:94 represent SKF and NASA versions of the code for three 
different shaft-bearing systems. The first set of examples 
prerent the results obtained with the SKF and NASA versions for 
a system in which an input pinion is supported by a preloaded 
pair of taper'^d roller bearings in a straddle configuration. 

The • 1 sd set of examples display program output for a 

system .. wnich an i. i-.\t pinion gear load is supported by a 
flanged cylindricf. Her beariug in conjunction with two 
angular contact ba . I bearings. 

The third set of examples represent execution for a single 
ball bearing;, sy:‘^en; operating under a combined radial and 
thrust load. 


TliO differences between the SKF and NASA versions of the code 
which are described in detail in Aopendix F, are reflected in 
the following output parameters: 

• Value of EHD film thickness. 

• Ratio of film thickness to surface roughness (H/Sigma) printed 
for most heavily loaded rolling element. 

• Lube-life factor. 

• Lio fatigue life. 

• Frictional heat generation rates. 

• Ratio of load carried by asperities to total contact load 
(Qasr/Qtot^ [for NASA version this ratio = 0, see 
Appendix F]. 
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NOMENCLATURE 
APPENDIX E 


16.9706n^U . 

^ 3 + 2k 

16.9706rioUjjR^ 

2 2 k 

(A + B )’ 

cage land radial clearance (in.) 
riolUyl K1 

cage pocket radial clearance (in.) 
R 



roller total length (in.) 
ball or roller radius (in.) 
cage land radius (in.) 



entrainment velocity (in. /sec.) 



original page is 


NOMENCLATURE continued 
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U 

i 


U 


u. 


V 


X 


V 


y 

Y 


Ho 


= entrainment velocity of the ring-land interface 
(in. /sec. ) 

= (o)^ + w^) for inner ring riding cage 

= Rj^ ((Do + for outer ring riding cage 

_ (in. /sec.) 

2 

_ V (in. /sec.) 

2 

~ (^x ^ ~ sliding velocity (in. /sec.) 

= sliding velocity at the ring-land interface (in. /sec.) 
= R£ (u)^ - (D^) for inner ring riding cage 
= R£ ((Do - (D^) for outer ring riding cage 
= R(Dy (in. /sec.) 

= -R(d (in. /sec.) 

- ' sin-^ (+) 

= angle between the rolling speed vector and the 
pumping force vector 

= absolute ambient viscosity ( lb. -sec. /in J) 
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yTHW. 

Oi ^ 

NOMENCLATURE continued 

APPENDIX 

(j) = cage orbital speed (rad. /sec.) 

c 

= inner ring orbital speed (rad. /sec.) 

Wo = outer ring orbital speed (rad. /sec.) 

= ball or roller rotational speed, x component (rad. /sec.) 
Wy = ball rotational speed, y component (rad. /sec.) 
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Calculation of Roller-to-Cage Pocket Forces 

A numerically stable means for calculating cage pocket forces 
in cylindrical roller bearings was developed in [24] . The analysis 
is generalized here to incluie cage pocket and cage land forces in 
ball and tapered roller bearings. 

Web geometry is taken as a radially outward cylindrical cavity 
for ball bearing simulations and a rectangular cavity for roller 
bearings, Figure £.1. The force exerted by a ball or roller is 
considered to act at the web midpoint, and all pocket contacts lie 
on the pitch circle. Interactions between roller ends and cage 
pocket sides are neglected. The load supported by the web is con- 
sidered to have two components; one normal to the plane of con- 
tact (the surface of the web) , and a frictional component in the 
plane of contact. 

The normal force component is calculated as a function of rol- 
ling element to cage pocket offset Z = Z ' where 

(jj. is the orbital speed of the i-th foiling element. Pocket loads 
ate computed as a function of the hydrodynamic lubricant film, 
h=h(Z ) that fills the gap between the cage web and the rolling 
elemeSt surface. 

Numerical stability is gained by taking a linear approximation 
to the force-displacement equation describing the roller-to-pocket 
normal load component 


^normal 


(E.l) 


where K is chosen such that equation (E.l) will match the exact 
solution at = 67N (15 lb.) for roller bearings, and 

^normal = IIN (2.51b.) for ball bearings. 


K = U 11 

“tc - .33Rno«^u)^ C 

^ Ir H 


(ball bearings) 


(E.2) 


K = 67 67 

(Cp - .33Rnoil^ ^ Cp 


(roller bearings) 
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Contact Inlet Region Hydrodynamic Friction Forces 
Ball /Cage Contacts 

In a properly lubricated ball bearing, a constant viscosity 
fluid film separates the metal surfaces of the ball and cage web. 
Surface motion causes pressure to build up at the entry to the 
contact. Ball surface motion entrains the lubricant into a 
gently narrowing wedge, creating a distributed traction on the 
ball surface. Area Integration of the surface traction yields a 
three dimensional load vector. Components of the load vector are 
given in terms of the dimensionless quantities = and ^ [25] . 

S 

Forces ; 

Rolling Components 

= %CoFj^ COSY (E.3) 

Fj^ = >jCoFj^(sinY) '\fk (E.4) 


Sliding Components 

^Sy ” (E.5) 

"sx = ^s''»''y‘Vy"» 

Moments ; 

‘ <^Ry ^ ^Sy' 

MFy = R • (Fj^ + Fg^) (E.8) 
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Values for F_ are obtained as follows. Prom [25] , the force 
directed normal co the plane of contact is 

''normal ‘ C.Pr^cosy 
or 


F_ = ^normal ^ (E.9) 

CoRjjCosy 

is known (E.l), and other terms are functions of geometry. 

norma 1 

Therefore, F„ is known. Functions shown in Figures E.2 and E.3 were 
curve fit* to give 

= = .26 F„ + 10.90 (E.IO) 

Fg R 


Roller/Caqe Contacts 

The analytic description of lubricated, rolle,*-to-cage pocket 
contact, Figure E.l (b) , is based on the lubrication of a rigid 
cylinder near a plane [26] . The radially directed friction force 
is expressed as a function of the normal force 


F 


f ric 


jj. F 

h norm 


(E.13) 


where is computed from (E.l). is a hydrodynamic friction 

coefficient 


* Fj^ is shown vs. dimensionless meniscus distance (pi) in Figure 
E.2; Fg is shown vs. Pi in Fig. E.3. The functions are closely 
approximated by: 


= 34.74 (tnpi) - 27.60 (E.ll) 


Fg = 8.82 (Jtnpi) + 3.89 (E.12) 

for hydrodynamic contact. We get (E.IO) by solving E.ll) for 
fn.^j and inserting this expression into (E.12). 
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norm 


Zn-U.a 


(E.14) 


The moment generated about the roller's x-axis is given by 


^x " ^^h ^norm 


(E.15) 


Cage Land Normal Forces and Friction Moment 

Forces which develop between a cage rail and its supporting 
ring surface are obtained using the hydrodynamic solution for 
self-acting, short journal bearings [27], Forces are assumed to 
act in the plane of rotation , 


The resultant of the pressure distribution on the cage can be 
described by orthogonal force components along and perpendicular to 
a line passing through the cage center and point of closest approach 
to the land. 


Figure E.4 shows tht relevant parameters for an inner ring land 

riding case, and Figure E.5 illustrates the outer ring land riding 

case.* The cage undergoes a radial displacement in the bearing 

XYZ frame, of magnitude e and direction 6^. An xyz frame is at- 

c 

tached to the cage, such that the y axis passes through the point 
of minimum film thicxness. 


Assuming an isoviscous, Newtonian fluid, the lubricant forces 
which develop at the guide ring are given as a function of ec- 
centricity [27] 



— P — 

n 


e 2 


(E.16) 


f; = ^ n»UgL- 


ne 


4(l-cM 


(E.17) 


* These figures were repeated from section 2.4.1. for convenience. 
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The upper sign applies to an inner ring and the lower to an ourer 
ring land riding cage. Eccentricity e is related to translation e, 

^ “ e (E.18) 

C 


For level 1 solutions » the radial force is set equal to the 

cage weight. The eccentricity £ necessary to support that weight 
is used to determine F^ and the cage-land friction torque 


m' « ± n V,RjL 2 IT 


(E.19) 


The angle is used to transform forces from the 
c 


frame to the cage reference frame: 



local xyz 


(E.20) 
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INTRODUCTION 


This appendix describes the SKF and NASA versions of 
SHABERTH. The primary differences between the two versions en- 
compass the EHD film thickness and the concentrated contact trac- 
tion force calculations. The relevant mathematical models are 
discussed in Section ”.2. The differences with respect to pro- 
gram input data are explained in Section F.3. 

F.2 MATHEMATICAL MODELS 
F.2.1 EHD Film Thickness 

In calculating the elastohydrodynamic film thickness, 
SHABERTH/SKF uses the Archard-Cowking equation (3) for point contact 
and the Dowson-Higginson equation (4) for line contact. Two film 
thickness reduction factors are then multiplicatively applied: 

1) a thermal factor due to heating in the 
contact inlet using the formulation of 
Cheng (11) 

2) a factor accounting for starvation at 
the contact developed by Chiu (12). 

SHABERTH/NASA uses the film thickness equation developed by 
Loewenthal et. al. (6). This equation is applicable for both 
point and line contacts. 

F . 2 . 2 Concentrat ed Con ta ct Tract ion 

The concentrated contact traction model used in SHABERTH/ 


F:2 



SKF accounts for lubricant shear and asperity interaction. 

A semi-empirical model developed by Chiu» discussed in (25), 
is used to calculate an EHD lubricant shear coefficient. 
Asperity effects are introduced by determining the portion 
of the contact load carried by the asperities, using the 
analysis of Tallian (5), and then calculating the resulting 
traction as the product of the normal load carried by the 
asperities times the asperity friction coefficient. In 
equation form the traction force is: 


^ ' ^EHD ^EHD ^ASP ’ '^ASP (1) 

^ ^ ^EHD ^ASP 


where 

F is the tracfion force 

^EHD is the normal load carried by the EHD film 

Ur-,,r> is the friction coefficient which develops 
EHD 

frojti lubricant shear 


'^ASP is the normal load carried by the asperities 
^ASr asperity friction coefficient 

Q is the total load 

SHABERTH/NASA calculates concentrated contact traction across 
the EHD film only, using the model developed by Allen, et. al. 
(7). This model determine* the traction force by first calcu* 
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lating the shear stress and then integrating the shear 
stress over the respective contact area. For a Newtonian fluid 
the shear stress is given by the equation 


T = n 




where 

X is the shear stress 
n is the dynamic viscosity 
V is the surface relative sliding velocity 
h is the film thickness 

The lubricant viscosity is assumed to be an exponential func- 
tion of pressure of the form: 


n = n^e 


a ’ s 


(4) 


where 


is the dynamic viscosity at atmospheric pressure 


rt is the pressure viscosity coefficient 
s is the normal stress 

The Allen formulation requires that the shear stress not exceed 
a specified fraction of the normal stress such that 


T = n 


V 

"fT 


if i 


T = Fs ifr — ■' Fs and nr > x 
h " ' 
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.1 


where 


T IS the critical shear stress for which a value of 
c 

2 

0.0069 N/mm (1000 psi) is normally used. 

T is called the lubricant friction coefficient and 
has been determined for specific lubricants. 

Typical values of T lie in the range 0.05 ^ F ^ 0.08. 

F.5 PROGRAM USE 

The selection of the desired .SHABTiRTH version has been 
made possible by the inclusion of two separate Map statements 
for the Univac 1100 computer: 

1) The original Map statement 

L3MAP.S AI.WAYS/MAP, ALWAYS/ABS 
enables execution of SHABERTH/SKF 

21 A neAv Map statement 

0MAP.S NASA/MAP, NASA/ABS 
enables execution of S11ABF,RTH/NASA 

The only difference between the two versions with respect 
to input data is on card type B16. For the NASA version, two 
additional lubricant data items are specified for NCODE ^ 0: 

1) AKN - Empirical lubricant constaht - columns 
71-75 

2) FRIG - Lubricant Friction Coefficient - columns 
76-80 
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Default values of AKN«5 0. 0 and FRIC-0.07 will be used if these 


spaces are left blank. For NCODE values of 1 to 4, the following 
values of AKN and FRIC are assigned. 


NCODE 

AKN 

FRIC 

1 

18.2 

0.075 

2 

18.2 

0.045 

3 

24.9 

0.070 

4 

18.2 

0.070 
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